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SUMMARY 

A 

Digital  synthetic  aperture  radar  (SAR)  imaging  techniques  have  pre¬ 
viously  only  been  reported  in  the  literature  in  a  fragmentary  manner.  This 
article  presents  a  comprehensive  review  of  the  theory  of  digital  SAR 
imaging  from  Earth  orbiting  satellites.  The  digital  SAR  imaging  process  is 
explained,  including  a  discussion  of  various  aspects  which  are  specific  to 
satellite-borne  SAR.  A  number  of  relevant  digital  processing  techniques 
are  reviewed  and  it  is  shown  how  these  techniques  may  be  applied  to  the 
processing  of  digital  SAR  data.  The  range  migration  problem  is  discussed 
and  various  techniques  for  overcoming  it  are  presented.  The  paper  should 
be  useful  not  only  to  the  designer  of  SAR  processors,  but  also  to  the  user 
of  digital  SAR  data,  and  images. 
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1  INTRODUCTION 

Synthetic  aperture  radar  imaging  of  the  Earth's  surface  was  proved  to  be  a 
practical  concept  by  the  SEASAT  global  ocean  monitoring  satellite  launched  by 
NASA  in  June  1978.  Although  this  satellite  remained  operational  only  until 
October  1978  a  considerable  quantity  of  SAR  digital  data  was  recorded,  and  is 
still  being  processed  at  various  centres  throughout  the  world. 

SEASAT  stimulated  great  interest  in  orbital  SAR  and  this  interest  is  grow¬ 
ing.  Since  1978  the  Shuttle  imaging  radar  SIR-A  has  provided  more  data,  although 
this  was  optically  recorded.  The  next  Shuttle  radar  mission  SIR-B  in  1984  will  be 
digitally  recorded,  and  looking  further  into  the  future  the  ESA  satellite  ERS-1, 
to  be  launched  in  late  1987,  will  also  carry  a  digital  SAR. 

Much  has  been  written  about  the  principles  of  SAR  and  of  optical  SAR  data 

1  2 

processing,  see  for  example  the  books  by  Harger  and  by  Hovanessian  and  the 
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collection  of  papers  edited  by  Kovaly  .  Further  background  information  can  be 
found  in  Refs  13-15.  The  literature  on  digital  SAR  data  processing  is,  however, 
somewhat  fragmented;  see  for  example  Refs  4-12.  The  SAR  processor  designer  needs 
to  have  a  theory  presented  as  a  coherent  whole.  Others  who  need  such  a  theory 
are  the  image  users  who  must  interpret  images  in  the  light  of  such  knowledge. 

The  material  presented  in  this  paper  is  in  the  nature  of  an  advanced  treat¬ 
ment  of  the  theory  underlying  orbital  synthetic  aperture  radars  and  provides  an 
introduction  to  the  techniques  of  digital  SAR  processing.  An  attempt  has  been 
made  to  present  the  theory  in  a  useful  form  as  a  coherent  whole..  In  addition, 
a  number  of  approximations  are  examined  and  an  attempt  made  to  answer  many  ques¬ 
tions  which  the  author  is  often  asked  by  workers  new  to  the  subject. 

An  introduction  to  the  basic  ideas  behind  SAR  will  be  found  in  Refs  1  and  2 
and  in  Refs  13-15.  It  is  assumed  in  the  treatment  presented  here  that  the 
reader  is  familiar  with  the  contents  of  these  papers. 

2  OUTLINE  OF  THE  SAR  IMAGING  PROCESS 

In  order  to  define  processes  and  terms  and  to  fix  ideas  an  outline  of  the 
SAR  imaging  process  is  presented  in  this  section.  In  particular,  as  an  illustra¬ 
tion,  the  image  of  a  single  point  on  an  absorbing  background  is  considered. 

Consider  just  one  pulse,  the  nth  pulse,  say,  and  an  echo  of  unit  amplitude 
from  a  point  on  an  absorbing  background.  The  received  echo  is: 

E(t)  =  exp  i2irff  (t  -  t  )  +  a(t  -  t  )2]  (1) 

u  u  n  D  n  D  J 


(see  Ref  10) . 
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If  the  radar  receiver  outputs  a  real  valued  signal  then,  ignoring  a 
constant  phase  term 

E(t)  =  cos  2ir[f  (t  -  t  )  +  a (t  -  t  )2j  ,  (2) 

fg  is  the  centre  frequency  of  the  transmitted  range  chirp,  t^  is  the  time 
measured  from  the  centre  of  the  nth  pulse,  t^  is  the  round  trip  time  from  radar 
to  point  to  radar  for  this  particular  pulse.  If  a(t)  is  the  range  of  the  point 
from  the  radar  then  t^  =  2a(t)/c  where  c  is  the  velocity  of  light.  In  this 
section  a(t)  is  assumed  constant  for  one  pulse  but  varies  from  pulse  to  pulse. 
This  is  the  'stop-start'  approximation  whereby  the  radar  platform  is  modelled  as 
being  stationary  while  a  pulse  is  being  transmitted  and  received;  t  is  the  over¬ 
all  time  coordinate. 

In  equations  (1)  and  (2)  the  transmitted  chirp  has  been  put  implicitly  into 
time  symmetric  form: 

A(t)  =  txp  i2irt  (f„  +  at  )  -T/2  <  t  <  T/2  (3) 

n  0  n 


-T/2  <  t  <  T/2 
n 


where  T  is  the  pulse  duration  and  a  =  bandwidth/2T  . 

A  number  of  interlinked  approximations  arc  implicitly  made  in  this  section. 
Briefly  these  are: 

(1)  The  'stop-start'  approximation  already  described. 

(2)  When  a  signal  is  received  as  a  result  of  scattering  from  a  target  and  there 
is  a  relative  velocity  between  transmitter/receiver  and  the  target  the  frequency 
of  the  received  signal  will  be  changed  by  a  scaling  factor  B  (the  Doppler 
effect) .  The  scaling  factor  B  is  taken  to  be  6=1+  2v/c  where  v  is  the 
relative  velocitv  of  the  radar  with  respect  to  the  target.  The  signal  spectrum 
is  transformed  s(f)  ■+  s(f/B)/»^|  6 1  and  f/B  is  taken  to  be  f  ( 1  —  2v/c)  and 
A&\  =  i  .  This  is  justified  by  the  fact  that  v  ^  c  . 

(3)  Although  for  a  spacecraft  SAR  the  signals  received  by  the  radar  were  trans¬ 
mitted  a  few  milliseconds  earlier  when  the  radar  was  at  another  position  it  is 
assumed  that  the  round  trip  range  and  Doppler  frequency  shift  are  those  corres¬ 
ponding  to  the  time  of  reception.  An  examination  of  this  approximation  for  the 
SEASAT  SAR  shows  that  the  main  effects  are  an  error  of  order  1  m  ir.  the  range 
direction  resulting  in  a  totally  negligible  defocussing  in  the  azimuth  direction 
and  an  image  shift  of  order  10  m  in  the  azimuth  direction  with  no  corresponding 
defocussing. 


These  approximations  are  associated  with  the  fact  that  during  the  time 
interval  for  which  a  point  on  the  ground  is  illuminated  by  a  pulse  (33.9  ys  for 
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SEASAT)  the  phase  change  caused  by  the  varying  range  is  very  small.  It  can  then 
be  shown  that  each  pulse  samples  the  phase  of  the  point,  that  is  4ira(t)/A  for 
the  nth  pulse.  This  whole  question  is  examined  in  section  3  where  it  is  shown 
that  the  'stop-start'  approximation  is  entirely  acceptable. 

2. 1  Range  compression 

On  receipt  of  each  pulse  the  radar  receiver  coherently  mixes  the  if  centre 
frequency  down  to  the  'offset  video  frequency'  f  .  The  offset  frequency  is 
usually  chosen  to  be  one  quarter  of  the  sampling  frequency  of  the  analogue-to- 
digital  converter.  The  reason  for  this  will  become  clear  in  section  6.1.  For 
t’EASAT  this  frequency  was  approximately  11.38  MHz.  The  received  echo  is  then: 

E(t)  =  cos  2tr[f  j  tn  -  fQtD  +  a(tn  -  tD)2]  .  (4) 

It  should  be  made  clear  at  this  stage  that  it  has  been  assumed  here  that  the 
signal  received  and  input  to  the  subsequent  processing  is  real-valued.  It  is 
entirely  possible  for  a  SAR  receiver  to  produce  a  complex  signd  with  a  single 
sideband  via  a  quadrature  filter  (see  Ref  16,  p  119).  In  that  case  an  offset 
video  frequency  of  zero  would  be  chosen.  However,  this  is  riot  usual  since  it 
is  easier  to  sample  the  real  valued  signal  and  convert  to  a  complex  valued  signal 
with  a  single  sided  spectrum  in  the  subsequent  digital  processing  -  see  section 
5.1. 


Range  compression  is  performed  on  each  pulse  by  correlating  the  pulse 
against  a  replica  of  the  transmitted  pulse  translated  tc  the  offset  video  band. 
The  peak  of  the  correlation  occurs  at  the  round  trip  delay  time  for  this  parti¬ 
cular  pulse,  ie  at  x  =  t  where  x  is  the  correlation  variable.  In  the 

K  JJ  K 

correlation  which  follows  the  time  symmetric  forms  of  transmitted  and  received 

pulses  are  used  and  the  start  of  transmission  occurs  at  t  =  -T/2  so  that 

n 

t^  =  0  coincides  with  the  centre  of  the  transmitted  pulse.  Correlating  (4) 
against  (3)  over  all  values  of  t^  from  t^  -  T/2  to  t^  +  T/2  gives  for  the 
range  point  spread  function 

tjj+T/2 

gR  =  *  /  ^  i2*[Vn  “  f0CD  +  a(tn  "  CD)2]  * 

tjj-T/2 

x  exp  -  i2^[fj(tn  -  tr)  +  a(tn  -  xR)2]dtn 
tD+T/2 

+  i  j  exp  -  i2*[f,tn  -  fQtD  +  a(tn  -  tj,)  ]  x 
tD-T/2 

exp  i2Tr[f,(tn  -  ?R)  ♦  o(tn  -  xR)2]dtn  , 


x 


(5) 
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and  so 

SR 


cos  2 irf-  fQtD  +  f  j  tr  -  «(tD 


sin  nAFR(tp  -  tr) 
"  TR) 


(6) 


where  AFR  =  the  chirp  (video)  bandwidth  =  2aT. 

Equation  (6)  assumes  continuous  functions.  In  reality  these  functions  are  sampled, 
and  the  correlation  is  performed  in  the  frequency  domain  via  an  FFT.  This  is 
explained  in  section  5.4.  The  result  ol  the  correlation  (6)  is,  of  course,  the 
same  if  the  correlation  is  performed  in  the  frequency  domain  via  Fourier  trans¬ 
forms,  although  -he  range  point  spread  function  (6)  is  slightly  different  when 
performed  on  sampled  data  and  is  periodically  repeated.  At  this  stage  of  the 
processing  the  cosine  term  in  equation  (6)  is  turned  into  an  exponential.  The 
reason  for  this  is  explained  in  section  5.1.  If  (6)  is  Fourier  transformed  with 
respect  to  v  one  obtains  an  approximately  rectangular  spectrum  (it  is  a  Fresnal 

integral)  of  width  AF  centred  on  f  multiolied  by  expi2nt  if  -f.n)  and  an 

R  I  i)  I  u 

approximately  rectangular  spectrum  of  width  AFR  centred  on  -f j  multiplied  by 
exp  -  i2irtD(fj  -fg)  •  As  was  mentioned  previously  the  range  correlation  is 
carried  out  in  the  frequency  domain  and  the  simple  step  of  setting  all  the  nega¬ 
tive  frequency  samples  to  zero  gives  only  the  positive  half  of  the  spectrum.  If 
this  is  inverse  Fourier  transformed  it  will  be  seen  that  instead  of  (6)  we  have: 


% 


exp  i2v[-  fQtD  +  fjtR  -  a(tD 


sin  7TAFR(tD  -  tr) 
^R^D  "  TR^ 


(7) 


At  the  same  time  the  offset  video  frequency  f j  is  mixed  to  zero  by  the  simple 
procedure  of  moving  all  the  frequency  samples  (see  section  8) .  The  resulting 
range  point  spread  function  is  then: 


g 


R 


exp  -  i2Tr[f0tD  +  a(tD 


sin  TTAFR(tp  -  tr) 
*AFR(tD  -  XR} 


<8) 


In  passing  it  should  be  pointed  out  that  these  expressions  for  g  have  been 

_1  K 

multiplied  by  a  constant  (=  A  z,  where  A  is  the  time  bandwidth  product)  in 
order  to  normalize  them. 


2.2  Azimuth  compression 


The  next  main  step  in  the  processing  is  azimuth  correlation  (although  in 
practice,  between  range  correlation  and  azimuth  correlation  there  are  the  very 
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important  steps  of  corner  turning  and  range  migration  correction.  These  are  com¬ 
putational  problems  and  are  examined  in  sections  2.3  and  10).  Azimuth  correla¬ 
tion  is  essentially  the  same  process  as  range  correlation  but  in  this  case  the 
frequency  modulation  results  from  the  Doppler  effect  and  the  change  in  t^  from 
pulse  to  pulse.  The  change  in  round  trip  time  t^  gives  a  rotating  vector  via 
the  phase  term  in  equation  (8) .  The  rate  at  which  it  rotates  represents  the 
Doppler  frequency.  In  a  sense,  then,  each  pulse  samples  the  phase  of  the  point 
and  over  a  number  of  pulses  the  rate  of  change  of  phase  is  built  up,  and  one  can 
then  meaningfully  refer  to  a  frequency.  The  change  in  phase  from  pulse  to  pulse 
is  a  linear  function  of  t^  in  equation  (8)  because  as  t^  varies  from  pulse  to 
pulse  the  value  of  x_  for  which  the  samples  are  selected  from  each  pulse  for 

IV- 

azimuth  compression  tracks  t^  with  an  arbitrary  offset.  The  quadratic  function 
of  tQ  in  (8)  therefore  remains  constant.  This  is  a  consequence  of  range  migra¬ 
tion  correction  and  is  explained  in  more  detail  in  section  2.3.  The  remaining 
phase  term  in  (8)  gives  the  pnase  history  of  the  point.  Now  we  may  write 


.  _  2a(t) 

D  c 


-!(■ 


+  *1'  *  v 


where  the  range  to  the  target  point  has  been  expanded  in  a  Taylor  series  about 
t  =  0  (ie  the  centre  of  the  synthetic  aperture).  That  is,  at  t  =  0  we  have 

a  =  a^  +  ajt  +  +  •••  OO) 

where  a^  is  the  slant  range  from  the  radar  to  the  point,  aj  is  the  slant  range 
velocity  and  a^  is  one  half  of  the  slant  range  acceleration.  Usually  only 
terms  up  to  the  quadratic  term  are  taken  into  account  but  sometimes,  for  extra 
precision,  the  cubic  term  is  used  as  in  the  RAE  processor.  The  range  polynomial 
(10)  is  calculated  for  the  general  case  of  an  elliptic  orbit  and  rotating  Earth 
in  section  4  where  it  is  shown  that  the  finite  value  of  aj  at  the  centre  of  the 
aperture  :s  a  result  of  radar  squint,  orbit  eccentricity  and  Earth  rotation. 


We  now  correlate  in  the  azimuth  direction  and  at  a  temporal  slant  range 
displacement  from  t  of  t  .  Note  that  t  is  not  constant.  This  is  the  so 
called  range  migration  effect.  Firstly,  because  the  Earth  rotates,  a  point  on  it 
is  taken  through  many  range  samples  (range  gates)  possibly  several  dozen  for  a 
high  resolution  spacecraft  SAR.  Secondly,  for  a  high  resolution  spacecraft  SAR 
the  geometry  is  such  that  an  arc  of  constant  slant  range  cannot  usually  be  con¬ 
sidered  'straight'  and  the  curvature  can  cause  the  point  to  move  through  a  number 
of  range  gates.  These  two  effects  are  known  as  range  walk  and  range  curvature. 
Together  they  constitute  range  migration  and  cause  serious  problems  for  the  SAR 
processor  designer.  These  effects  are  considered  further  in  sections  2.3  and  13. 
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non 


In  passing,  note  that  the  range  time  coordinate  t  is  not  constant,  it 

tracks  tQ  and  has  an  arbitrary  constant  offset  from  it  corresponding  to  the 

displacement  of  the  line  of  data  being  used  for  azimuth  compression  from  the 

centre  of  the  range  point  spread  function  | t  -  t  |  . 

D  R 

Some  approximations  are  now  made  in  order  to  make  the  azimuth  correlation 
process  clearer  although  it  should  be  noted  that  these  are  not  actually  made  in 
a  working  digital  processor.  For  a  small  change  in  slant  range  the  first  order 
changes  to  the  range  polynomial  ate  Aa^  in  a^  ,  Aa^^/a^)  in  a^  and 
Aa^^/ a^)  in  a ^  .  For  lengths  comparable  with  the  point  spread  function 
scale  it  can  easily  be  shown  that  the  changes  to  aj  and  a^  are  negligible, 
so  that  only  the  change  Aa^  in  a^  ,  need  be  considered.  Thus  expanding  a 
about  a  =  a^  and  t  =  0  the  following  terms  are  retained: 

a'(t)  =  aQ  +  AaQ  +  ajt  +  a2t2  (11) 

where  Aa^  =  tr  -j  .  (12) 

In  addition,  for  a  displacement  d  in  the  direction  of  satellite  motion,  there 

will  be  a  shift  in  the  azimuth  time  coordinate  x.  equal  to  d/v  where  v 

A  p  p 

is  the  local  velocity  of  the  satellite. 

Here  is  ample  scope  for  confusion  here.  Conventionally,  the  azimuth 

direction  is  defined  as  the  direction  of  the  satellite  motion.  This  has  the 

advantage  of  being  constant  over  the  whole  swath  width.  Strictly,  however,  it 

ought  to  be  defined  as  the  direction  of  the  resultant  of  v  and  da(t)/dt  at 

P 

the  relevant  point.  This  latter  definition  is  used  here  and  the  two  directions 
are  called  the  'along  track  direction'  and  the  'azimuth  direction'.  The  true 
azimuth  direction,  then,  varies  over  the  swath  and  if  an  image  of  a  bright  point 
is  observed  it  will  be  seen  that  the  sidelobes  in  the  azimuth  direction  are 
aligned  along  a  slight  curve  which  does  not  intersect  the  range  direction  at 
right  angles.  This  is  illustrated  in  Fig  1 . 

The  actual  phase  history  of  the  point  is  exp(-  ^wf^t^)  with  t^  given 
by  (9)  and  the  phase  history  against  which  it  is  correlated  by  the  processor  is 
exp(-  i2Trfgt^)  with  t^  given  by: 

CD  =  U30  +  Aa0  +  3l(t  +  V  +  a2(t  +  TA)2J  *  (,3) 
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It  should  be  noted  that  these  various  approximations  have  been  introduced  here  in 

order  to  make  this  explanation  clear.  In  pract:ce  digital  processors  compute  the 

phase  history  directly  from  the  spacecraft  orbit,  the  Earth  rotation  and  the 

geometry.  Correlating  (8)  against  exp(-  i2irf  t')  ,  then,  gives: 

U 


8<tr’ta> 


exp  -  -  tr) 

T./2 


2  Sin  tAFR< fcp  ~  V 
lrAFR(tD  '  V 


x  /  eXp  -  1  -T  2f0[a0  +  V  +  a2t2]  X 


-T./2 
A 


2f0  r 

X  exp  i 2tt  laQ  +  AaQ  +  a^t  -  xA>  +  a2(t 


TA)2]dt 


<1A) 


where  a  synthetic  aperture  interval  of  ±T^/2  centred  on  zero  has  been  taken. 
An  asymmetric  aperture  simply  introduces  an  additional  phase  term.  Sc 


2f 

2  0  /  2\ 

g(xR,TA)  =  exp  -  i2ira(tD  -  xR)  exp  12tt  fAa0  +  a^  +  a^J  x 


sin  ^AFR(tD  -  xR)  sin 

irAFR(  fcD  ~  V  WAFATA 


(15) 


AFa  is  the  Doppler  bandwidth.  The  point  spread  function  8(tr»t^^  can 

expressed  in  lengths  in  azimuth  (x)  and  slant  range  (y)  by  means  of  x.  =  x/v 

A  p 

and  x  =  2y/c  .  The  point  spread  function  in  (15)  has  been  scaled  by  multiply- 

K  _i  _i 

ing  by  A  2B  where  A  and  B  are  the  range  and  azimuth  time  bandwidth 
products. 

2.3  Range  migration 

It  will  be  observed  that  the  position  of  the  range  point  spread  function 
gR  (see  equation  (8))  is  a  function  of  time  because  t^  is  a  function  of  time 
as  it  must  be  to  obtain  a  Doppler  frequency  shift  on  which  the  synthetic  aperture 
principle  depends.  If  this  shift  is  much  smaller  than  the  range  sample  separation 
then  its  effect  on  the  azimuth  correlation  integral  can  be  ignored:  for  example 
most  aircraft  synthetic  aperture  radar  systems  are  designed  so  that  this  is  the 
case.  However,  as  pointed  out  in  the  last  section,  it  is  not  possible  to  do  this 
in  the  case  of  an  orbital  SAR  for  a  high  resolution  imaging  system.  The  effect  of 
the  temporal  dependence  of  the  range  point  spread  function  is  to  couple  the  range  and 


along  track  directions  in  the  azimuth  correlation  integral.  If  the  azimuth 
correlation  integral  is  performed  along  the  'true'  azimuth  direction  then  the 
range  and  azimuth  correlations  are  uncoupled.  The  'true*  azimuth  direction  was 
defined  in  the  last  section  in  one  way,  and  in  this  section  it  is  implicitly 
defined  in  another,  but  closely  related  and  equivalent  way.  Returning  to  the 
azimuth  correlation  (14)  and  putting  the  time  dependent  range  point  spread  func¬ 
tion  inside  the  integral  we  obtain 


e^R’XA* 


fA  . Sin  "4FR(tD  -  TR> 

j  6XP  -  l2"“(tD  -  V  .AF-(tn  -  -  )  - 

rf  /  o  K.  JJ  K 


c  0 1_ 


aQ  +  ajt  +  a2t  x 


u  r  2 

x  exp  i2ir  — -  |aQ  +  AaQ  +  a,(t  -  ta)  +  a2(t  -  ta> 


Jdt  . 


The  range  and  along  track  components  are  uncoupled  by  performing  an  additional 
convolution  operation  along  the  line  defined  by 


xo  * !( 


°0  *  V  +  V 


in  range  sample  time  and  azimuth  time  space.  Thus  if  t_  =  t  -  t  : 

K  D  K 


8(VtA>  =  f 


T./2  T  +T/2 
f  A  D  f 


6(tD  ‘  CR  -  TR)dTR  * 


-T./2  t  -T/2 
A  D 


sin  ttAF  (t  -  t  ) 

X  exp  -  i2„0,(tD  -  tr)  — fR7t— —  « 


2ir  T 

exp  -  i  —  2fQ jaQ  +  a,t 


+  V  * 


X  exp  i 2 7t  jaQ  *  AaQ  +  a,(t  -  ta>  +  a2(t  -  xA)2]dt 


where  the  range  compressed  data  spans  the  time  interval  tD  ±  T/2  .  It  will  be 
noted  that  the  additional  operation  must  be  carried  out  before  azimuth  correlation. 


1 1 

Usually  the  convolution  is  carried  out  by  'data  selection'  using  more  or  less 
approximate  methods;  this  is  explained  further  in  section  14. 

3  THE  STOP-START  APPROXIMATION 

In  section  2  the  range  compression  process  and  subsequent  azimuth  com¬ 
pression  were  based  on  the  'stop-start'  approximation  in  which  the  delay  time 
tp  was  taken  to  be  constant  over  a  pulse  width.  This  approximation,  and  others 
closely  linked  with  it,  outlined  in  section  2,  are  now  examined  in  more  detail  in 
this  section.  In  digital  processing  of  any  kind  and  in  particular  in  digital  SAR 
processing  it  is  necessary  to  have  a  very  clear  analysis  of  all  the  details  cf 
the  imaging  process  including  the  approximations  since  a  mathematical  model  must 
be  constructed.  Failure  to  do  this  can  lead  to  unsatisfactory  results  and 
unexpected  problems. 

3. 1  The  delay  time 

Let  the  radar  be  at  rest  in  a  reference  frame  S  with  coordinates 
(x,y,z,t)  and  let  the  target  point  be  at  rest  in  a  frame  S'  with  coordinates 
(x' ,y' ,z' ,t') .  There  is  a  relative  velocity  between  the  frames  and  the  two 
frames  are  related  via  a  Lorentz  transformation  (see  Ref  17,  p  1 1 C  and  also 
p  132).  For  the  moment  any  accelerations  are  ignored.  A  radar  operates  by 
sending  out  a  pulse  of  electromagnetic  waves  at  time  tj  say.,  and  receiving  the 
pulse  at  t^  both  times  being  measured  in  frame  S  at  rest  relative  to  the 
radar.  The  time  at  which  the  pulse  is  estimated  to  arrive  at  the  target  point 
by  an  observer  at  the  radar  is  then  (t^  +  tj)/2  ;  hence 

t2  =  (t,  +  t3)/2  .  (18) 

A  different  time  is  measured  by  the  observer  at  the  target  since  to  him  it 

appears  that  the  clocks  in  the  radar  reference  frame  are  slow.  In  the  Lorentz 

transformation  which  connects  an  event  in  each  frame  the  relative  velocity 

2  2. 

appears  as  v  /c  .  An  Earth-orbiting  spacecraft  may  have  a  velocity  of  several 

kilometres  per  second,  say  6.7  km/s  at  a  (typical)  orbital  height  of  800  km  and 
2  2  -10 

v  / c  —  5  x  10  .  This  is  generally  not  significant  for  a  low  frequency  SAR 

such  as  the  L-band  SEASAT  but  note  that  at  a  wavelength  of  3  cm  the  carrier 
frequency  is  10*^  Hz  and  hence  the  second  order  term  in  v/c  leads  to  a  Doppler 
shift  of  5  Hz  or  so  which  might  be  measurable.  Thus  with  future  X-band  or 
higher  frequency  orbiting  SARs  the  effect  may  have  to  be  taken  into  account  as 
will  be  explained  later  in  the  section. 

2  2 

In  the  analysis  which  follows  v/c  is  very  small  and  so  the  factor 
/l  -  v^/‘.2  in  the  Lorentz  transformation  from  one  frame  to  another  will 
always  be  taken  as  unity.  Then  the  analysis  will  be  correct  up  to  first 
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order  terms  in  v/c  ,  Within  this  approximation,  time  can  be  taken  as  absolute 
and  the  analysis  can  be  performed  in  'absolute'  space. 

Let  the  position  vector  of  the  radar  be  r^ (t)  and  that  of  the  target 
point  ?2(t)  in  any  suitable  coordinate  system.  An  electromagnetic  wave  leaves 
the  radar  at  time  tj  and  travels  towards  the  target  irriving  there  at  time  t 2  > 
is  reflected  and  arrives  back  at  the  radar  at  t^  .  The  distance  travelled  by 
the  wave  from  radar  to  target 

dj  =  |?j (t)  -  ?2(t) j  (19) 

and  from  target  to  radar 

d2  =  |r2(t)  -  rj(t)|  .  (20) 

Hence  the  total  time  of  the  wavefront  from  radar  to  target  to  radar  is 

c3  -  tj  ~  ^|r2(t2)  ~  ri^t3^l  +  “  *2^2^]  ' 

Since  we  are  dealing  with  differences  it  is  easier  to  redefine  the  coordinate 
system  so  that  the  origin  is  at  the  target  point,  ie  ?0(t)  s  0  .  Then 

?j(t)  =  r (t) 

and  so 

t3  -  tj  =  -^  [r(t2)  +  r(tj)]  .  (22) 

In  equation  (22)  the  time  of  reception  corresponding  to  t^  is  t  and  the 
delay  time  is  tp  .  Hence 

CD  -  7  ['«*„>  +  <23> 


tp  must  now  be  obtained  explicitly  by  solving  this  equation.  Fortunately  t^  is 
usually  small  (around  6  ms  for  a  slant  range  of  850  km)  and  ?(t)  is  a  smooth 
function  so  that  r(t^  -  t^)  can  be  expanded  as  a  Taylor  series  around  r(rn)  • 
Now, 


r  (t 

n 


tD) 


f(tn)  -  v(tn)tD  +  Ktn)t2/2 


(24) 


where  v  =  3r/3t  and  I  =  3v/9t  ,  which  gives 


r(tn  -  V 


r(t  )  -  r(c  ' 
n  n 


<*(4) 


r  being  the  unit  vector  in  the  d  rection  of  r  .  It  is  not  necessary  to 
calculate  the  second  order  tel*m  since  the  eventual  expression  for  t^  wil 


which  gives 


■  i  [2r(tn>  -  i(\)  ■  * 0 (•‘dJj 


y>y 

!  3 

■S&. 

(25) 

■  /  >_ 

-  -  v 

only 

(25) 

<< 

(26) 

U 
»  *.*  *. 

\v;: 

(27) 


A  more  elaborate  analysis  including  the  second  order  terms  confirms  this  equation 

2  2 

and  shows  that  the  acceleration  terms  are  of  order  v  /c  and  thus  negligible. 

The  angle  between  r(t)  ana  v(t)  is  the  azimuth  angle  <J>  ,  equal  to  ir/2  when 

the  radar  is  broadside  on  to  the  target  point.  One  might  expect  the  term 

r(t  )  •  v(t  )/c  to  be  negligibly  small  and  this  turns  out  to  be  the  case, 
n  n 

We  have 


2r(t  )  [*  v(t  )  /  2v1 

T2-  [»  '  IT"  COS  +  °(t)J  * 


(28) 


Note  that  at  this  stage  the  time  delay  tn  is  still  a  function  of  time  t^  via 


r(t  ),  v(t  )  and  <J»(t  )  . 
n  n  n 


Define  to  be  the  time  at  the  centre  of  the  received  pulse  when  t^  is 
assumed  constant,  ie  if  v  *  0  and  corresponds  to  the  fixed  tQ  of  the  stop-start 
approximation.  The  delay  time  at  Tq  is  then 


2r( 


W 


v  r 

=  I  2c 


cos  ^>(Tq)  +  0 


Q1  • 


(29) 


We  require  the  time  delay  at  time  offsets  from  up  to  ±T/2  .  The  time  delay 

at  time  t  is  then 


tD(T)  *  tD(T0)  +  tD(T0)T  +  +  ' 


(30) 


Differentiating  (28)  to  give  the  coefficients  in  this  Taylor  expansion  about  t 

r 

gives  terms  involving  2r(tn)/c  and  its  derivatives  and  vr(t^)  cos  <f>(tn)/c 
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and  its  derivatives.  The  maximum  rate,  of  change  of  <j>  occurs  for  the  broadside 
mode  and  is  typically  a  few  tens  of  mrad/s  maximum.  Also,  as  always,  v  ^  c  and 
the  leading  terms  only  are  taken  in  the  expansion  to  give 


2r(T0)  -  2?<t0) 
VV  =  and  tD(Tc)  =  — ~ 


and  hence 


■*[■ 


tD(T)  =  c  r(T0)  +  r(T0)T  + 


<Sir 


;>,v| 

-VS 


3 . 2  Range  compression  with  varying  time  delay 


.  -i.-i 
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where 


tD(T)  "  T0  +  T1T  +  Xel 


-fr +  0(T } 
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where 


and  from  (31) 


In  addition 


tD(x)  =  X0  +  2VlX  +  Xe2 


2X0X2X  .  2  2  .  n,  3s 

- -  +  t  t  +  0(x  )  , 


2r(x0) 


2r(V 


2r (xq) 


c  *  '2 


h  =  j  exp  ’  “"[Vd  ‘ 

-T/2 


f 1 XR  +  a(tD  '  XR)  +  2aT(tD  -  XR} 


•Jdr  , 


This  is  obtained  by  rearranging  the  first  integral  in  (5),  the  second  is  ignored 
for  the  purposes  of  the  present  analysis.  Substituting  from  (32)  and  (34)  we 
obtain 
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bR 


=  exp  -  i27r[f0x0  -  f,xR  -  «(tr  -  xQ)  2]  x 


T/2 

x  J  exp  -  i2irjf0 


x  .  -  ax  0  +  2ax.x  +  2axx  .  +  2ax_x  .  x 
el  e2  1  el  0  el 


-T/2 


,] 


f*  f  T  " 

0  1 

X  exp  -  i2ir  2ax  (xQ  -  xR)  + 


dx 


(37) 


The  first  exponential  under  the  integral  contains  quadratic  and  higher  powers  of 
x  and  the  second  exponential  contains  the  linear  terms  in  x  .  The  non-linear 
powers  of  x  result  in  defocussing  and  the  linear  terms  simply  result  in  a  shift 
in  the  position  of  the  point  spread  function  in  the  range  direction.  The  quadra¬ 
tic  phase  term  is  slowly  varying  since  the  non-linear  terms  are  small.  The 
linear  phase  term  is  also  slowly  varying  in  the  neighbourhood  of  the  peak  of  the 


point  spread  function  since 


T  —  T 

0  R 


near  the  peak.  The  non-linear  term  is 


f  x  i- 
0  2  2 


+  ax.x2(2  -  x.)  +  0(x3) 


(38) 


Note  that  x,  =  2r(xQ)/c  2  since  r  ^  v  and  v  ^  c  .  Also  2aXj  ^ 

since  a  is  typically  of  the  order  of  10 12  Hz/s  .  The  largest  non-linear 

2 

term  is  then  2aTjX  .  In  section  4  it  is  demonstrated  that  a  typical  value 
of  r  for  a  broadside  looking  orbital  SAR  at  an  altitude  of  800  km  is  of 
order  100  ro/s.  The  quadratic  phase  change  over  the  correlation  is  then  of 
order  1  mrad  for  a  pulse  width  of  the  order  of  tens  of  microseconds.  This 
magnitude  of  phase  change  is  completely  negligible. 

The  first  order  efrect  of  the  stop-start  approximation  is  a  shift  in  the 
position  of  the  point  spread  function  by  fgXj/2a  ,  which,  for  SEASAT,  is  about 
1.5  ns  or  0.45  m  in  slant  range.  This  error,  however,  could  be  much  greater  for 
a  radar  operating  at  a  higher  frequency.  Also,  in  a  SAR  operating  in  a  squint 
mode  t j  could  easily  be  an  order  of  magnitude  greater.  A  shift  in  the  posi¬ 
tion  of  the  point  spread  function,  however,  causes  no  real  problems  unless  the 
shift  varies  over  the  synthetic  aperture  by  more  than,  say  one  half  resolution 
cell,  tj  does  vary  in  the  azimuth  direction  but  not  by  enough  to  cause 
problems  with  the  approximation  considered  here.  Thus  in  conclusion  the 
stop-start  approximation  is  valid  to  a  high  degree  of  precision  for  orbital 
synthetic  aperture  radars. 
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4  THE  RANGE  POLYNOMIAL 

In  section  2.2,  eouation  (10),  the  range  from  radar  to  point  was  expanded 
as  a  function  of  time: 

2  3 

a(t)  =  aQ  +  ajt  +  a2t  +  a.jt  ...  .  (10) 

In  this  section  the  terms  in  the  expansion  are  calculated  for  a  general  elliptic 
orbit,  and  various  approximations  considered. 

It  seems  appropriate  to  mention  here  that  there  are  at  least  three  modes  in 
which  a  SAR  may  operate.  These  are  the  normal  mode,  the  squint  mode  and  the 
spotlight  mode.  In  the  normal  mode  the  radar  looks  out  sideways  to  the  track  and 
perpendicular  to  it.  The  squint  mode  is  similar  except  that  the  radar  is  perma¬ 
nently  squinted  either  forwards  or  backwards  to  the  normal  mode.  In  the  spot¬ 
light  mode  the  radar  squint  angle  is  varied  continuously  as  the  radar  passes  a 
target  so  that  the  target  is  constantly  illuminated,  and  hence  the  synthetic 
aperture  length  is  no  longer  limited  by  the  footprint  width.  A  very  large  aper¬ 
ture  can  then,  in  principle,  be  synthesised.  The  squint  angle  is  varied  by  vary¬ 
ing  the  pitch  angle  and/or  the  yaw  angle  of  the  spacecraft  depending  on  the  angle 
of  incidence  of  the  radar. 

The  analysis  which  follows  is  applicable  (within  limits)  to  all  three 
modes.  Before  calculating  the  range  polynomial,  however,  some  considerations 
concerning  the  choice  of  orbit  for  SAR  spacecraft  are  presented. 

4. 1  Orbits  for  SAR  satellites 

Circular  or  nearly  circular  orbits  are  arguably  the  best  choice  for  radar 
remote  sensing  satellites.  There  are  several  reasons  for  this.  First,  if  an 
approximately  constant  height  is  maintained  above  the  Earth's  surface  then  the 
Doppler  bandwidth  will  remain  constant  and  this  leads  to  a  constant  pulse 
recurrence  frequency.  The  prf  is  one  of  the  main  design  parameters  for  a 
synthetic  aperture  radar  and  if  it  can  be  kept  constant  this  greatly  simplifies 
the  design  of  the  radar  system.  Further  discussion  of  this  point  belongs  to  the 
design  of  the  radar  rather  than  the  processing  and  is  not  pursued  here. 

Secondly,  if  the  radar  data  are  required  to  be  processed  into  images  soon 
after  the  satellite  pass  then  nominal  (and  possibly  inaccurate)  orbital  elements 
will  have  to  be  used  in  the  processing.  Refined  orbital  parameters  are  not  avail¬ 
able  until  several  days  have  elapsed  because  of  the  necessity  of  measuring  the 
orbit  at  several  ground  stations  around  the  Earth.  The  inaccuracy  in  the  pro¬ 
cessing  due  to  the  use  of  nominal  orbital  elements  seems  to  be  minimised  if 
a  circular  orbit  is  employed,  or  if  the  orbit  is  frozen  in  some  way. 


It  could  be  argued  that  it  might  be  desirable  to  have  a  slightly  elliptical 
orbit  so  as  to  maintain  (to  first,  order)  a  constant  height  above  the  Earth.  The 
orbit  apogee  would  then  lie  on  the  equator,  the  Earth  being  a  flattened  spheroid 
bulging  at  the  equator.  This  argument  was  advanced  in  Ref  19  in  connection  with 
the  SEASAT  altimeter  and  such  an  orbit  has  an  advantage  because  if  the  argument 
of  perigee  is  90°  (which,  of  course,  it  is  if  the  apogee  lies  on  the  equator) 
then,  in  principle,  it  is  possible  to  choose  an  inclination  and  eccentricity 
such  that  the  orbit  is  frozen.  That  is  to  say  such  that  the  perigee  does  not 
precess.  Normally,  the  perigee  processes  in  the  orbit  plane  by  a  few  degrees  per 
day  due  to  the  Earth's  oblateness.  The  reader  is  referred  to  Ref  19  for  a  fuller 
account . 

Another  effect  of  the  Earth's  oblateness  is  to  precess  the  plane  of  the 
orbit  by  a  few  degrees  per  day.  The  rate  of  precession  ft  is  given  approxi¬ 
mately  by:  ?^2 

ft  =  -  9.97^— )  — " -°-S"j  ~2  degrees/day  (39) 

3  (1  -  e  ) 

see  Ref  18,  equation  1.1.  It  will  be  observed  that  this  is  a  function  of  semi¬ 
major  axis  a  ,  eccentricity  e  and  angle  of  inclination  i  . 

This  effect  has  some  importance  to  spacecraft  designs  because  the  rate  of 

precession  can  be  chosen  to  be  such  that  the  solar  panels  which  power  most 

satellites  continuously  face  the  Sun.  This  is  a  Sun-synchronous  orbit  and  has 

a  rate  of  precession  of  approximately  360°  in  365  days.  Not  too  much  weight 

should  be  placed  on  this  requirement,  however,  because  quite  a  large  variation 

on  a  perfect  Sun-synchronous  orbit  is  acceptable  and  not  all  spacecraft  are 

powered  by  solar  cells.  If  one  were  designing  for  a  synchronous  orbit  then  once 

the  height  of  the  satellite  above  the  ground  is  chosen  (this  is  a  function  of 

radar  maximum  power,  area  of  coverage,  maximum  data  rate  capacity  etc)  then  the 

inclination  angle  i  follows  for  a  given  eccentricity.  The  required  inclination 
•  o  .  • 

is  greater  than  90  since  0  must  be  positive  and  such  orbits  are  known  as 
retrograde  orbits.  For  example  in  the  case  of  SEASAT  i  was  slightly  greater 
than  108°.  Of  course,  if  the  radar  platform  is  not  dependent  on  the  Sun  for  its 
power  (as  is  the  case  far  the  Shuttle  for  example)  then  the  orbit  is  not  con¬ 
strained  in  this  manner. 

The  above  discussion  is  a  somewhat  simplified  account  of  the  choice  of 
orbit  for  a  SAR  satellite.  The  actual  details  are  complicated  and  further 
details  of  the  choice  of  orbit  for  SEASAT  will  be  found  in  Ref  19. 
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In  the  analysis  which  follows  a  typical  orbit  is  assumed  to  be  as  in  Fig  2. 
The  orbit  is  shown  projected  on  to  a  spherical  surface  which  is  ’fixed'  relative 
to  the  background  stars.  The  direction  labelled  T  is  the  vernal  equinox  (the 
intersection  of  the  Earth's  equator  and  the  plane  of  the  ecliptic)  and  the  centre 
of  the  surface  coincides  with  the  centre  of  the  Earth.  The  point  P  is  the  point 
on  the  Earth's  surface  for  which  the  range  polynomial  is  to  be  computed.  The 
latitude  of  the  point  is  measured  in  geocentric  coordinates  and  is  constant,  but 
the  longitude  is  constantly  increasing  due  to  the  Earth's  rotation.  The  point  A 
is  the  projection  of  the  real  antenna  focus  on  to  the  celestial  sphere. 


4.2  Latitude  and  longitude  changes  of  radar  and  target 

Let  the  point  P  in  Fig  2  have  a  geocentric  latitude  and  longitude  measured 
on  the  fixed  surface  of  i b  and  <p  .  Let  R  be  the  distance  of  the  focus  of  the 

(real)  radar  antenna  from  the  Earth's  centre  and  r  the  distance  of  P  from  the 

Earth's  centre.  Then  if  a  is  the  distance  between  the  radar  antenna  focus  and 
the  point  P 

a2  =  R2  +  r2  -  2Rr  cos  6  (40) 

where  0  is  the  angle  AOP  =  AOD  in  Fig  2.  0  is  a  function  of  time  due  to  the 

motion  of  the  satellite  and  the  rotation  of  the  Earth  as  indicated  in  Fig  2. 

R  is  also  a  function  of  time  due  to  the  non-circular  o_uic,  but  r  is  constant. 


In  Fig  2  D  is  the  projection  of  r  on  to  the  celestial  sphere  upon  which 
the  orbit  has  been  projected,  and  the  arc  AD  passing  through  the  projection  of 
the  antenna  focus  A  and  D  has  on  a  great  circle  centre  0  .  The  meridians 
passing  through  A  and  D  are  also  great  circles  and  triangle  ADN,  N  being 
the  north  pole  is  a  spherical  triangle.  Let  the  latitude  and  longitude  of  A  be 
¥  and  $  in  geocentric  coordinates.  Then  from  AADN  we  have: 


cos  0  =  sin  T  sin  ip  +  cos  ¥  cos  cos($  -  <j>) 


(41) 


which  gives  0  in  terms 

The  purpose  of  the 
in  the  range  polynomial 

a 


of  the  latitudes  and  longitudes  of  A(4',$)  and 


analysis  here  is  to  compute  the  coefficients  a 
2  3 

=  a0  +  a, t  +  a2t  +  a3t  ... 


o 


in  the  neighbourhood  of  t  =  0  .  This  involves  differentiating  (40)  successively 
and  setting  t  =  0  .  Hence  explicitly  we  have: 


»q  =  a  at  t  =  0 
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(42) 


aoai 


R0(Rq  -  r  cos  60)  -  R0r  ^  (cos  0Q) 


(43) 


2a0a2  +  al  =  R0^R0  "  r  COS  V  +  Ro[R0  ”  2r  dt  ^LS  6(P_ 


"  V  TT  (cos  60) 

at 


(44) 


6aQa3  +  6aja2  =  RQ(R0  “  r  cos 


V  +  RoK  "  3r  IE  (cos  60)] 


+  fto[So  - 


_£ 

dt: 


>] 


3r  — ^  (cos  eQ)|  -  RQr  — j(cos  0Q)  . 


(45) 


The  subscript  (0)  indicates  that  the  variable  is  taken  at  t  =  0  .  When  the 
orbit  is  circular  R^  and  higher  derivatives  are  zero  and  the  equations  simplify 
greatly.  To  proceed  further  it  is  necessary  to  compute  derivatives  of  cos  6 
from  (41)  in  the  region  of  t  =  0  . 


Let  Fig  2  show  the  positions  of  radar  and  point  P  at  t  =  0  ,  and  let 
the  projection  of  the  satellite  move  on  to  position  B  at  time  t  ,  see  Fig  3. 
In  Fig  3  C  is  a  point  on  the  celestial  sphere  which  has  the  same  geocentric 
latitude  as  B  and  the  same  longitude  as  A  (the  position  of  the  radar  antenna 

A 

focus  at  t  =  0).  The  angle  BAC  =  and  is  the  'local  heading  angle'  of  the 
satellite  at  t  =  0  .  The  angle  AOB  is  the  angle  through  which  the  satellite 
has  moved  in  its  orbit  from  A  to  B  and  AOB  =  A a  .  Also  BOC  =  AT  and 

A  A 

AOC  =  Af  .  The  change  in  longitude,  A$  =  BNA,  AT  and  A$  are  related  by 


.  AT 
sin  -j- 


.  A4> 

sin  -7 


cos^  +  A^) 


(46) 


The  analysis  is  somewhat  intricate  and  for  that  reason  a  number  of  the  inter¬ 
mediate  steps  are  now  set  down. 

The  arc  BC  lies  on  a  great  circle  and  so  AABC  is  a  spherical  triangle. 

A 

The  angle  ACB  is  not  exactly  a  right  angle,  although  provided  Aa  is  small 
(which  it  always  is  even  for  the  'spotlight'  modes)  then  ACB  is  very  nearly  a 
rignt  angle.  For  example  Aa  is  usually  of  the  order  of  1  mrad  for  normal  modes 
Angle  ACB  is  therefore  written  as  tt/2  +  e  .  From  the  sin  formula  for 
spherical  triangles 
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AT  =  sin 


in '(' 


sin  Vq  sin  Aa\ 
cos  e  , 


To  find  e  consider  ACBN  in  which  NCB  =  tt/2  -  e  and  then 


=  tan"1!"  sin(Y  +  AY)  ilZ™*  Ml]  . 

L  sin  A$  J 


Spherical  triangle  ABC  gives 


AY  =  tan  tan  Aa  cos  - 


sin  AF  sin  e 
cos  Aa  cos  AY 


which  is  an  implicit  formula  for  AY  . 

It  will  be  evident  that  computing  the  coefficients  in  the  range  polynomial 
directly  L<y  means  of  the  above  formulae  is  not  a  practicable  task.  Fortunately, 
however,  this  is  not  necessary.  All  the  angular  changes  Aip,  A$,  Aa,  AT  and  t 
are  cf  the  same  order  of  magnitude.  Also  Aa  is  a  first  order  function  of 
time.  All  terms  in  the  range  polynomial  up  to  and  including  the  third  order 
term  a^  are  required  and  hence  it  is  only  necessary  to  expand  the  expressions 
for  Ay,  e,  AY  and  up  to  third  order  to  obtain  the  required  polynomial 

coefficients  exactly.  Expanding  (46)  gives: 


A$  =  Ar  sec  Y  +  AYAr  sec  Y„  tan  Yn  +  AY2AT  tanV.  sec  Y. 

U  0  0  0  0 

AY2AF  ...  at3  2  3  3 

+  — sec  Yq  +  sec  YQ  tan  YQ  +  0(ATAY  ,  AYAr  )  (50) 

(47)  gives: 

3  2 

Ar  =  Aa  sin  vQ  -  sin  vQ  cosVq  +  ~~  sin  vQ  +  C(Aa5 ,  Aa3t2,Aae4)  (51) 


(48)  gives: 


o 

A$  .  AYA$  a*  A4 

- sin  +  ^  cos  *0  ’  sin  *0 


A$3  .2  33 

+  “24"  sm  Yq  cos  YQ  -  0(AYA$  ,A$AY  ) 


(52) 


V>  *i|,.«  w  1 


and  (49)  gives: 

3 

AY  =  Aa  cos  Vq  -  eAT  +  py  cos  Vg  sin^Vg  +  0(Aa^Are,A4'^Are,Ar^e,Ar^)  « 


(53) 


These  equations  can  be  manipulated  to  give,  finally: 

2  3 

....  .  Aa  .  2  ...  Aa  .  3 

AT  =  Aa  cos  Vq - y  sin  Vg  tan  Yg  +  -y—  cos  Vq  sin  Vg 

3 

-  yp  sin2Vg  cos  \)q  sec^Tg  +  O(Aa^)  (54) 

2 

A$  =  Aa  sin  Vq  sec  Yg  +  Aa  sin  Vg  cos  Vg  sec  Yg  tan  Yg 

Aa^  2  3  2  2 

3“  sin  Vg  cos  Vg  sec  Yg  +  Aa  sxn  Vg  cos  Vg  tan  “g  sec  Yg 

3 

-  yp  sin^Vg  tan^Yg  sec  Yg  +  O(Aa^)  ,  (55) 


which  are  the  required  formulae  for  the  changes  in  latitude  and  longitude  in 
going  from  A  to  B  in  Fig  3,  expressed  in  terns  of  change  in  satellite  angle 
(or  'true  anomaly’)  Aa  . 


A  problem  arises  here  due  to  the  singularity  at  the  poles  where  Yg  =  ±ir/2 
and  it  is  necessary  to  avoid  this  case  in  the  analysis.  Thus  the  region  within 
about  2  of  the  north  and  south  poles  is  excluded.  Fortunately  this  presents 
no  difficulty  because  an  exactly  polar  orbit  is  very  unlikely  to  be  used  in 
practice.  The  suffix  (0)  on  the  variables  in  equations  (54)  and  (55)  refer  to 
t  =  0  ,  ie  mint  A  in  Fig  3. 

Differentiating  (54)  and  (55)  and  setting  t  =  0  gives: 


Yg  =  a0  cos  Vg 


(56) 


*o  =  “o  cos  vj  "  “o  sin2vo  t3n 


(57) 


“o  cos  vo  "  ?V:J  sin2vo  tan  vo  +  2“o  cos  vo  sin% 


,•3  .  2  2W 

-  3dg  Sin  Vg  COS  Vg  sec  Yg 


(58) 


.V.^vH  A7*^.  ^Vr^t v.V.Vk 


a^  sin  Vq  sec  ¥q 


(59) 


*•  *2 

ciq  sin  Vq  sec  4^  +  2a^  sin  Vq  cos  Vq  sec  Yq  tan  ¥q 


(60) 


o0  sin  Vq  sec  4^  +  60^0^  sin  Vq  cos  sec  4^  tan  u'q 

+  2aQsinVp  cos^Vq  sec  ^  +  6a^  sin  cos^v^  tan^  sec  4^ 


-  2a^  sin^Vp  tan^^  sec  Vq 


(61) 


In  addition  to  the  motion  of  the  satellite  we  also  have  to  consider  the  rotation 
of  the  Earth.  The  latitude  of  a  given  point  P  remains  constant  but  its  longi¬ 
tude  increases  at  a  constant  rate  5  rad/s.  Thus 


-  4 


(62) 


At  this  point  the  idea  of  the  squint  angle  is  introduced.  The  squint  angle  is 
here  defined  as  the  angle  between  the  great  circle  arc  AD  and  the  great  circle 
perpendicular  to  the  orbit  direction  at  A  -  see  Fig  2.  This  definition  is 
valid  for  general  elliptic  orbits  as  well  as  circular  orbits  since  we  consider 
the  projection  on  to  a  sphere.  If  the  squint  angle  is  zero  at  the  centre  of  the 
synthetic  aperture  at  t  =  0  then  the  radar  is  operating  in  the  normal 
’broadside’  mode.  This  angle  could  be  (and  generally  is)  a  function  of  time  due 
to  charges  in  the  pitch  and  yaw  angles  of  the  spacecraft.  The  squint  angle  is 
defined  to  be  positive  when  looking  forwards. 

A 

Angle  NAD  (Fig  2)  is  then  tt/2  “  vq  “  °q  where  Oq  is  the  squint  angle  at 
t  =  0.  Note  that  throughout  this  section  is  defined  as  positive  in  Figs  2  and  3. 

We  are  now  in  a  position  to  compute  the  coefficients  in  the  range 
polynomials. 

4.3  The  linear  term 

Differentiating  (41),  substituting  (56)  and  (59)  and  simplifying  by  means 
of  two  spherical  triangle  formulae  for  triangle  ADN,  viz 

sin  cos  'Fq  cos(Vq  ►  Oq)  =  sin  4^  cos($q  -  <$q)  cos  4>q  cos(Vq  +  Oq) 


+  sin($Q  -  4>q) 


sin(vQ  +  Oq)  cos  (63) 
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and 


sin  0Q  cos(vQ  +  aQ)  =  cos  i|»0  sin(<J>Q  -  <j)Q) 


(64) 


gives 


d (cos  eQ) 

- -  =  aQ  sin  eQ  sin  aQ  +  £  sin  0Q  cos  4^  cos(vQ  +  aQ) 


(65) 


which  together  with  (43)  gives: 


al  =  ^  (R0  “  r  cos  eo) 


V  •  . 

—  aQ  sin  60  sin  aQ 


Ror 

—  5  sin  0Q  cos  ¥0  cos(vQ  +  aQ) 


(66) 


Thus  the  linear  (range  walk)  term  neatly  separates  into  three  terms.  The  first 
is  zero  for  a  circular  orbit  since  R  is  then  zero,  and  is  thus  associated  with 
orbit  eccentricity.  The  second  term  is  zero  when  the  squint  angle  Oq  is  zero 
and  is  thus  associated  with  the  radar  squint.  The  third  term  is  zero  when  £  is 
zero  and  is  thus  associated  with  the  rotation  of  the  Earth.  Thus  linear  range 
walk  is  in  general  caused  by  three  factors;  orbit  eccentricity,  squint  and  Earth 
rotation. 


It  is  instructive  to  calculate  a^  for  a  typical  example  and  the  example 
chosen  (Fig  4)  is  an  image  from  Seasat  Rev. 762  over  the  east  coast  of  England 
(the  Wash  and  East  Anglia) .  The  coefficient  aj  will  be  calculated  for  the 
centre  of  the  image.  The  Keplerian  orbital  elements  for  Seasat  Rev. 762  are  given 
in  Table  1 ; 

Table  1 

Orbital  elements  for  Seasat  Rev. 762 


Semi -major  axis 

7161 .39494  km 

Eccentricity 

0.00186004 

Inclination 

108.0203° 

Argument  of  perigee 

148.1649°  ' 

Right  ascension 

89.3670° 

Mean  anomaly 

252.3242° 

I 


£ 


r 


S 


vV 


V 

V* 

rj* 


$ 

a 


•r.  . 

% 


m 


I 


£ 


,v, 


9* 


^T 
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Rq  and  its  derivatives  and  cXq  and  higher  derivatives  can  easily  be  worked  out 
from  the  orbital  elements  using  the  usual  orbit  equations  -  see  for  example 
Ref  18,  Chapter  3. 

The  latitude  and  longitudes  of  the  satellite  nadir  and  the  centre  of  the 

•  ••  • 

image  at  t  =  0  (tne  centre  of  the  aperture)  and  values  of  r,  R^,  R^,  R^,  a^, 

(Xq  and  0q  are  given  in  Table  2. 

r  is  the  Earth  radius  at  P  and  is  calculated  on  the  basis  of  a  spheroidal 
Earth  with  a  flattening  factor  of  1/298.3  and  equatorial  radius  of  6.37816  y  lO^m. 


Table  2 

Parameters  of  centre  of  Fig  4  at  t  =  0 


Latitude  of  nadir 

51.25° 

Longitude  of  nadir 

-2.08° 

Latitude  of  point 

52.78° 

Longitude  of  point 

0.67° 

o 

CD 

2.265° 

r 

6.36444  x  106  m 

Ro 

7.16297  x  106  m 

Ro 

1 .37760  x  10  m/s 

Eo 

»o 

-1.758803  x  l(f3  m/s2 
1.041312  x  10'3  rad/s 

ao 

4.005  x  io-9  rad/s2 

5 

0.7292115  x  io”4  rad/s 

vo 

29.20° 

The  corresponding  values  of  the  three  terms  in  the  expression  for  are: 

*0 

Eccentric  term:  —  (Rrt  -  r  cos  dn)  =  13.147  m/s 

aQ  0  0 


Squint  term: 


Rotation  term: 


V  .  . 

- an  sm  6- 

aQ  0  0 


V 

-  5  sin  9_  cos  cos 

Sq  U  U  U 


-38.887  m/s  per  degree  of  squint 
for  small  squint  angles 

-82.385  m/s  . 


The  pitch  and  yaw  angles  of  the  satellite  of  the  centre  of  the  aperture  were 
about  -0.0040  rad  and  0.0112  rad.  Because  of  the  steep  (low)  incidence  angle 
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the  pitch  angle  has  a  greater  effect  on  the  squint  than  the  yaw  and  in  this 
instance  they  nearly  cancel  out  t'.  give  a  squint  angle  of  about  -0.092°.  Hence 
the  squint  term  amounts  to  some  3.58  m/s.  The  squint  angle  can  easily  be  0.5°  or 
so  on  occasions  for  SEASAT. 

The  linear  (range  walk)  term  is  then  about  65.66  m/s  for  this  example. 

One  can  draw  some  very  interesting  conclusions  from  equation  (66).  Notice 
that  it  is  possible  to  select  a  squint  angle  Oq  such  that  aj  is  zero,  aence 
linear  range  walk  could,  in  principle,  be  eliminated  for  a  specific  value  of  0q 
corresponding  to,  say,  the  middle  of  the  swath  by  varying  the  squint  angle  of 
the  radar.  Moreover  if  the  orbit  is  circular  then 


RqT  sin  0 


sin  Oq  r  -| 

T -  [“o  sin  au  +  ?  cos  Yo  cos(vo  +  Vj 


and  a^  would  be  zero  over  the  whole  swath  width. 

In  practice  the  centre  of  the  Doppler  spectrum  (which  is  amplitude 
modulated  by  the  antenna  beam  pattern)  would  be  measured  on  board  the  satellitf 
and  the  satellite  yaw  or  (better)  the  pitch  angle  would  be  varied  via  the 
satellite  attitude  control  system  until  the  Doppler  spectrum  is  centred  on  zero. 
This  type  of  technique  could  have  far  reaching  implications  for  the  design  of  a 
SAR  and  processor  system  for  real-time  processing.  The  range  walk  problem  is 
explained  further  in  section  12. 

4.4  The  quadratic  term 

Differentiating  (41)  twice,  substituting  from  (56),  (57),  (59)  and  (60)  and 
simplifying  by  means  of  (41)  and  an  additional  formula  for  triangle  ADN: 


sin  0q  sin(vQ  +  aQ)  =  sin  cos  YQ  -  cos  sin  YQ  cos($0  -  (f>Q)  (68) 


gives 


^  •  2  •• 

— 2  (cos  0q)  *  -  Oq  cos  0Q  +  Oq  sin  0Q  sin  aQ 

dt 


2£a0Jsinv0  cos  ipQ  cos($Q  -  <f>Q)  -  cos  vQ  cos  sinYg  sin(«Q  -  <|>0)J 


-  C  cos  Yq  cos  t0  cos($0-<f>0)  . 


The  terms  including  £  can  be  simplified  further  but  only  by  introducing  yet 
another  angle  and  so  it  has  been  left  as  it  is.  In  conjunction  with  (43)  and 
(45),  (69)  gives: 
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_l  r  +i!o 

2ao  0  Ro 


J  [e0  -  r  cos  e0]  +  2^  [r„  -  of]  - 


i(or  r  *2  .. 

['  °0  COS  0o  +  a0  Sin  90  sin  °0 


+  2^o{sin  v0  COS  ^0  C0S<VV  -  cos  vo  cos  ^0  s*n  sin($o" v} 


-  5  cos  Vq  cos  ipQ  cos($q  -  4>q) 


which  is  a  general  expression  for  the  quadratic  (range  curvature)  term  for 
elliptic  orbits  and  rotating  Earth.  It  can  be  separated  into  terms  associated 
with  eccentricity  (the  terms  containing  derivatives  of  R  and  second  derivative 
of  a),  and  those  associated  with  Earth  rotation  (the  terms  containing  £ ) .  For 
a  circular  orbit,  stationary  Earth  and  zero  squint  angle  the  expression  simpli¬ 
fies  drastically  to  give 


Ror“o  cos  9o 


2ao  V  V 


cos  0, 


where  is  the  satellite  velocity  =  RgCtg  and  h  is  the  height  of  the 

satellite  above  the  Earth.  Since  0^  is  usually  small  (2°  to  3°)  and 
h  —  Rq/10  then  the  usual  simple  expression  given  in  the  literature  for  a„, 
Vp^2a0  ’  See  ^°r  examPle  Re^  equation  (1),  is  correct  to  within  10%  or  so. 

It  .s  again  instructive  to  determine  the  magnitudes  of  the  various  terms  in  (70) 
for  the  'typical  example'  Fig  4.  Substituting  values  from  Table  2  into  the  part 
of  (70)  which  contains  the  radial  derivatives  gives: 


f  2R21 

_l  r  +  _2  r 

%  L 0  ao  J  t 


Rq  -  r  cos  8q 


]  .  1  r^2  2l  R0al 

J  *  ^  l*»  -  *ij  "C-  =  - 


3.795  x  10'3  m/s2  . 


The  remainder  of  (70)  contains  the  angular  derivatives: 


1ST  "  ao  cos  eo  +  ao  sin  eo  sin  ao 


+  25a0|sin  vQ  cos  ipQ  cos(0Q -4>Q)  -  cos  vQ  cos  i|>Q  sin  4'0  sin($0  -  <)>0)| 
-  5^  cos  4*^  cos  cos($q-^)q)^  =  28.0946  m/s^  .  1 


.0946  m/s 


If  the  maximum  synthetic  aperture  time  is  taken  at  ±1  s,  and  if  the  maximum 

tolerable  phase  error  over  the  aperture  is  taken  to  be  equivalent  to  A/ 16  where 

A  is  the  wavelength,  then  it  is  obvious  that  the  terms  containing  the  radial 

•  ♦  »  » 

derivatives  are  negligible.  The  term  containing  a.  is  also  negligible 

“5  2  ^ 

(R^ra^  sin  6q  sin  o^/2a^  =  4.86  x  10  m/s  ). 

There  is  no  doubt  that  the  terms  associated  with  eccentricity  can  be 
ignored  in  this  instance  and  that  a ^  can  be  calculated  on  the  basis  of  a 
circular  orbit.  Hence  such  phenomena  as  the  precession  of  the  perigee  and  varia¬ 
tion  of  eccentricity  have  no  material  effect  on  the  quadratic  term  -  provided 
that  the  orbit  eccentricity  is  small  enough.  A  more  refined  analysis  can  set 
limits  on  the  orbit  eccentricity  but  this  would  take  us  too  far  from  the  main 
theme  of  this  paper. 

4.5  The  cubic  term 

The  terms  involving  eccentricity  in  the  quadratic  case  were  shown  to  be 
negligible  and  they  are  even  more  so  in  the  cubic  case.  Unfortunately 
all  the  other  terms  (of  which  there  are  very  many)  in  the  equation  for  a^  only 
differ  by  one  order  of  magnitude  or  so  and  all  must  therefore  be  considered. 

The  resulting  expression  for  a^  will  not  be  set  down  here;  it  can  easily  be 
calculated  along  the  same  lines  as  for  aj  and  a ^  .  The  expressions  are 
complicated  and  not  easily  simplified  due  to  the  inherent  asymmetry.  However  if 
this  is  done  it  will  be  found  that  substituting  the  values  in  Table  2  for  the 
example  given  before  gives  a^  =  1.26  x  10  m/s  .  Hence  over  a  synthetic  aper¬ 
ture  of  ±1  s  there  is  a  variation  of  about  1/9  wavelength.  This,  unfortunately, 
is  on  the  limits  of  the  maximum  acceptable  error,  and  to  be  on  the  safe  side  in 
this  case  the  cubic  term  should  be  included.  The  example  is  not,  of  course,  a 
worst  case  example  and  the  cubic  term  increases  at  more  northerly  latitudes. 

For  SAR  systems  with  shorter  wavelengths  (such  as  the  proposed  ERS-1  SAR 
with  A  —  5.6  cm)  and  hence  shorter  synthetic  aperture  intervals,  the  cubic  term 
will  be  completely  negligible.  Conversely  with  systems  which  have  a  long 
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wavelength  and  long  synthetic  aperture  interval  (such  as  the  proposed  spotlight 
mode  for  the  Shuttle  SIR-B  mission)  the  cubic  term  will  be  significant. 

4.6  The  range  polynomial:  conclusion 

In  this  section  it  has  been  shown  how  it  is  possible  to  calculate  the 
terms  in  the  range  polynomial  for  the  most  general  case  of  an  elliptical  orbit 
and  rotating  Earth.  It  has  been  shown  that  the  eccentricity  of  the  radar  plat¬ 
form  orbit  has  most  effect  on  the  linear  term,  for  (small  eccentricity  it  can  be 
neglected  so  far  as  the  quadratic  and  cubic  terms  are  concerned.  A  more  refined 
analysis  is  possible  based  on  this  section  but  will  not  be  presented  here.  One 
indication  of  this  is  that  the  use  of  nominal  orbital  parameters  is  more  likely 
to  have  an  effect  on  the  linear  range  walk  correction  in  the  processing  than  on 
focussing  via  the  quadratic  phase  term.  The  use  of  ’auto-focussing'  techniques 
to  overcome  the  effects  of  the  use  of  nominal  orbital  elements  is  thus  somewhat 
questionable.  All  auto-focussing  techniques  known  to  the  author  rely  in  some  way 
on  the  defocussing  being  caused  by  an  error  on  the  quadratic  phase  over  the 
correlation.  An  error  on  the  linear  range  walk  can  also  cause  defocussing  due  to 
the  selection  of  incorrect  data  in  the  range  walk  correction  procedure  and  this 
is  discussed  again  in  section  13. 

5  DIGITAL  PROCESSING:  MATHEMATICAL  PRELIMINARIES 

In  this  section  some  relevant  background  theory  is  briefly  presented.  Most 
of  this  theory  will  be  found  in  one  form  or  another  in  Ref  16.  Consider  a  func¬ 
tion  of  time  f(t)  ;  f(t)  may  be  complex  valued.  Its  Fourier  transform  F(co) 
is 

00 

F(o>)  =  J 
—00 

The  inverse  transform  is 

F(t)  =  ~ 

Ztt 


f(t)  exp  *  iwt  dt  .  (74) 


00 

j  F(w)  exp  iuit  dus  .  (75) 

— 00 


F(w)  is  the  spectrum  of  the  function  f(t)  ,  and  is  a  complex  valued  function. 
If  f(t)  is  purely  real  or  imaginary  then  the  real  and  imaginary  parts  of  F(tu) 
have  symmetry  properties: 


if  f(t)  is  real  valued  then 
if  f(t)  is  imaginary  valued  then 


F*(w)  =  F(-w) 

(76) 

F*(to)  =  -  F(-oj) 

(77) 

n 


*  % 
"Cv 


•'s' 


i.* . 


where  the  asterisk  indicates  complex  conjugate. 
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5. 1  The  analytic  signal 

From  equations  (76)  and  (77)  we  see  that  a  real  (or  imaginary)  function 
can  therefore  be  completely  determined  by  only  the  positive  (or  negative) 
frequencies,  ie  only  half  of  the  spectrum  because  the  other  half  can  be  recon¬ 
structed  via  equatiors  (76)  and  (77).  This  takes  us  to  the  idea  of  the  analytic 

signal. 

Let  f(t)  be  purely  real  with  F(to)  given  by  (74).  Then  F*  (to)  =  F  (-to) . 
Let  g(t)  be  purely  imaginary  with  Fourier  transform  G(<d)  .  Then 
G*  (w)  =  — G  (  — ui)  . 

Let  g(t)  be  such  that  for  m  >  0  ,  G(o>)  =  F(m)  . 

Then  G*(u>)  =  F*(w)  as  well,  for  m  >  0  . 

Then  for  to  <  G  we  have  G(-oo)  =  -F(oi)  . 

Let  z(t)  =  f(t)  +  g(t)  ,  then  its  Fourier  transform  is 

F(u))  +  G(w)  =  2F(a>)  w  >  0 

=  0  a)  <  0 


and  we  therefore  have  a  single  sided  spectrum.  Usually  g(t)  is  written  as 
if(t)  and  then 

z(t)  =  f(t)  +  if(t) 


z(oj)  = 


2F(w) 

0 


0)  >  0 
a)  <  0 


And  so  starting  with  a  real  function  f(t)  we  can  form  its  Fourier  transform 
F(co)  ,  double  it,  and  set  its  negative  frequencies  equal  to  zero.  If  this 
function  is  then  inverse  transformed,  the  real  part  of  the  resulting  complex 
valued  function  z(t)  is  equal  to  the  original  real  function  f(t)  .  Thus  only 
half  of  the  spectrum  of  the  real  valued  f(t)  is  needed  to  define  it  completely. 


Now  consider  the  function 


a(t)  exp  i<p ( t)  =  a(t)  cos  (t)  +  ia(t)  sin  <Kt) 

where  a(t)  and  <j>(t)  are  real  valued.  a(t)  exp  i<Kt.)  is  assumed  to  be 
analytic  (and  is  an  entire  function  of  exponential  type  with  singularities  only 
at  infinity).  Then  it  can  be  proved  (see  Ref  20,  Chapter  5)  that  the  function 
has  a  single  sided  spectrum.  Hence  it  is  possible  to  identify  a(t)  cos  4>(t) 
with  f(t)  and  a(t)  sin  4> (t)  with  f(t)  .  In  actual  fact  f(t)  and  f(t) 


H  ^  tC^  * 


are  mutual  Hilbert  transforms,  but  that  is  only  of  passing  interest  here.  It 
will  be  observed  that  the  chirp  signal,  equation  (3),  is  in  this  form  and  is  an 
analytic  signal.  Likewise  (1)  is  the  analytic  signal  corresponding  to  the  real 
signal  (2). 

Thus  if  we  have  a  real  valued  function  such  as  a(t)  cos  <p( t)  writing  it 
as  a(t)  exp  iq> ( t)  gives  the  ’analytic  signal'  with  a  single  sided  spectrum. 

It  will  be  shown  in  section  6.1  that,  in  the  case  of  a  SAR  system  in  which 
the  Doppler  band  extends' over  both  positive  and  negative  frequencies,  it  is 
absolutely  necessary  to  operate  on  analytic  signals. 

5.2  The  convolution  and  correlation  theorems 


g(t)  =  /  f  (r)h(t  -  T)dx  =  /  f(t  -  i)h(-r)dT 


G(u))  =  F(u))H(m) 


(see  Ref  16,  p.21,  equation  (1-47)).  The  convolution  is  usually  written 


g(t)  =  f  (t)  ©  h(t)  . 

Hence  convolution  in  the  time  domain  corresponds  to  multiplication  in  the 
frequency  domain  (and  vice  versa) . 

A  similar  theorem  exists  for  the  correlation  function: 


f*(r)h(t  +  T)di 


f*(r  -  t)h(r)dr 


C(u)  -  F*(aj)H(to) 


In  Ref  21,  equation  (406)  p.768,  it  is  shown  that  the  complex  response  from  a 
matched  filter  is  determined  by  a  correlation  process. 

To  perform  matched  filtering  in  the  frequency  domain  therefore,  one 
Fourier  transforms  the  signal  and  multiplies  it  by  the  complex  conjugate  of  the 
chirp  replica  spectrum  and  then  performs  the  inverse  transformation.  This  is  the 
basis  of  frequency  domain  compression  in  both  the  range  and  azimuth  directions. 


The  sampled  discrete  Fourier  transform 

The  discrete  Fourier  transform  pair  analogous  to  the  continuous  transforms 


are: 


FU) 


N-l 


f(k)  exp  - 


i2Trk£ 

N 


k=0 
N-l 


f(k) 


£=0 


F(£)  exp 


i2TrkH 

N 


(83) 


(84) 


This  pair  of  transforms  can  be  regarded  as  a  mapping  of  the  sequence  f(k), 
k  =  0  ...  N-l  on  to  the  sequence  F(£),  1  =  0  ...  N-l  and  vice  versa.  There 
are  N  samples  in  the  'time'  domain  and  N  in  the  'frequency'  domain. 

The  reader  is  reminded  that  it  is  implicit  in  this  representation  that  the 

original  continuous  time  function  f(t)  and  the  corresponding  spectrum  F(w) 

have  both  been  periodically  repeated.  Usually  the  temporal  function  f(t)  has 

a  finite  length  and  is  sampled  at  intervals  of  At  ,  say,  so  that  if  T  is  the 

total  length  T  =  NAt  .  The  corresponding  spectrum,  however,  is  not  finite  and 

the  process  of  periodically  repeating  it  causes  the  edges  to  be  periodically 

folded  into  the  fundamental  period.  This  is  called  aliasing.  The  width  of  the 

fundamental  frequency  period  is  1/At  ,  see  Ref  22,  and  it  is  therefore  necessary 

that  there  shall  be  no  frequencies  higher  than  f  /2,  where  f  =  1 /At,  present 

s  s 

in  the  original  spectrum.  The  reason  for  the  additional  factor  of  J  is  that  we 

have  both  the  positive  and  negative  frequencies  to  contend  with  (see  Fig  5). 

The  frequency  fg/2  is  commonly  called  the  folding  frequency.  It  will  be 

observed  that  the  sample  sequence  F(&)  runs  from  0  to  N  -  1  and  therefore  in  the 

first  half  of  the  sampled  spectrum  we  have  positive  frequencies  and  in  the  second 

half  negative  frequencies.  The  spectrum  is  folded  about  f  /2  . 

s 

There  are  no  aliasing  problems  with  the  time  sampled  sequence  because  it  is 

a  finite  length  record  and  is  zero  outside  the  sampled  range.  Each  sample 

in  the  time  domain  is  separated  by  At  where  At  =  T/N  and  each  sample  in  the 

frequency  domain  is  separated  by  Af  =  1/T  where  Af  =  F/N  and  F  =  f  the 

s 

sampling  frequency. 

This  discussion  applies  to  a  real  valued  function.  When  the  sampled  func¬ 
tion  is  complex  valued  with  a  single  sided  spectrum  then  the  spectrum  can  be 
periodically  repeated  at  double  the  rate  of  a  double  sided  spectrum  without  ambi¬ 
guity.  This  means  that  if  the  sampling  frequency  is  f  the  spectrum  can  have 

s 

frequencies  up  to  f  instead  of  f  /2  before  ambiguity  sets  in.  This  can  be 

s  s 
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illustrated  by  a  simple  example.  Consider  the  real  and  complex  signals  cos  2nft 


and  exp  i2irft  .  Suppose  both  the  signals  are  sampled  at  times  n/fs,  n  =  0  ...  N  -  I, 


Then  for  the  complex  signal  frequencies  f  and  f  +  f  are  ambiguous  because 


exp 


i2u(f  +  f  )n 
s 


i2irf 


=  exp 


n 


(85) 


and  for  the  real  signal  the  frequencies  f  and  f  +  fs/2  are  (just)  ambiguous 


because 


cos 


r,  <f+V«l 

I  2tt  - j -  =  -  cos  — 

L  s  J  s 


(86) 


If  real  signals  with  frequencies  in  the  range  0  to  f  /2  are  sampled  at  fn 


then  the  complex  signal  may  be  sampled  at  fg/2  without  ambiguity. 


5.4  The  discrete  convolution  and  correlation  theorems 


If 


N-l 


gOO 


N-l 


f  (m)h(k  -m) 


m=0 


■il 


f  (k  -m)h(m) 


(87) 


m=0 


then 


G(Jl)  =  FU)HU) 


l  =  0  ...  N  -  1.  (88) 


The  convolutions  in  equation  (87)  and  (88)  are  cyclic  because  f(m)  and  h(m) 
m  =  0  ...  N-l  are  periodic  functions.  Again,  then,  convolution  in  the  sampled 
time  space  corresponds  to  multiplication  in  the  sampled  spectrum  space. 


There  is  a  corresponding  result  for  the  correlation. 


N-l 


If 


c(k) 


N-l 


f*(m)h(k  +m) 


m=0 


f*  (m  -  k)h(m) 


(89) 


m=0 


then 


c  (A)  =  F*(£)H(t)  , 


(90) 


again,  this  correlation  is  cyclic. 
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Cyclic  convolutions  can  produce  problems  (see  Ref  22,  p.48).  For  example 
in  the  case  of  the  range  compression  process  it  will  be  demonstrated  that  there 


are  end  effects  which  make  it  necessary  to  process  more  data  than  is  required  by 
the  image. 

5.5  Interpolation  via  the  Fourier  transform 


In  terms  of  its  discrete  Fourier  transform  (DFT)  a  sampled  function  is: 


N-l 

f(kAt)  =  ^^F(2.Au>)  exp  ~ 

1=0 


0  <  k  <  N  -  1 
0  <  £  <  N  -  1 , 


(see  equation  (84)).  At  is  the  spacing  of  the  time  samples 
spacing  of  the  angular  frequency  samples.  Am  =  2irAf  '  here 
sample  spacing 


1  _  _2tt_ 

NAf  NAw  * 


and 

Af 


Aw  is  the 
is  the  frequency 

(92) 


Suppose  that  the  length  of  the  inverse  DFT  is  doubled  and  define  a  new  set  of 
spectrum  samples  F'(fcAw)  such  that 


F'(AAu)  =  F(£Au>) 
=  0 


for  0  <  l  <  N  -  1 
for  N  <  l  <  2N  -  1 . 


(93) 


Then  the  new  inverse  DFT  is: 


2N-1 


(£Aw) 


i2ir£k 
exp  — 2N~ 


F(£Aw)  exp 


£=0 


i2ir£k 

2N 


N-l 

x  i2ir£(k/2) 

=  >  F(£Au)  exp  - - - 

£=0 


(94) 


Comparing  this  with  equation  (91)  gives: 

2N-1 

^  F'  (*Aw)  exp  =  f  ~  .  (95) 

£=0 


Since  the  number  of  samples  in  the  inverse  DFT  was  do-'bled,  the  numbers  of 
samples  of  the  time  function  f(t)  was  also  doubled  and  there  are  therefore 


'.•.’A 
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2N  samples  at  intervals  of  At/2  .  That  is  to  say  the  time  function  has  been 
interpolated.  All  the  original  samples  exist  and  in  addition  there  are  new 
(interpolated)  samples  midway  between  the  original  samples. 

This  process  can  evidently  be  continued  and  if  the  inverse  DFT  is  per¬ 


formed  with  a  length  2 
points. 


times  the  forward  DFT  this  gives 


2n  -  1 


in  erpolation 


Some  further  comments  will  be  made  on  the  process  in  section  6.3.  There  is 
another  interpolation  technique  using  the  DFT  which  is  not  limited  to  dividing 
the  sample  interval  by  some  power  of  2.  Returning  to  equation  (91)  it  is  evident 
that 

N-l 

f(kAf  +  aAt)  =  ^F(£Aw)  exp  exp  .  (96) 

£=0 


It  is  therefore  possible  to  produce  new  samples  at  any  chosen  point  within  the 
original  intervals  (a  <  1)  simply  by  weighting  the  spectrum  samples  by  a  phase 
shift  exp(i2ira£/N)  before  inverse  Fourier  transforming. 

Both  of  the  above  techniques  have  successfully  been  used  on  SaR  processing 
at  RAE. 

6  SPECTRA 

6. 1  The  range  spectrum 

The  signal  transmitted  by  a  SAR  is  usually  a  frequency  modulated  'chirp' 
with  a  large  time-bandwidth  product.  The  spectrum  of  such  a  chirp  is  calculated 
in  Ref  16  (example  8-5,  p.270).  It  is  a  Fresnel  integral.  If  the  chirp  is 
centred  on  f^  (the  radar  'carrier'  or  centre  frequency)  and  has  a  nominal 
width  of  AF  (=  the  chirp  rate  in  Hz/s  *  the  pulse  width)  then  the  spectrum  is 
as  sketched  in  Fig  6.  This  is  the  spectrum  of  the  real  valued  signal  (2)  as 
obtained  by  the  radar  receiver  from  a  point  scatterer.  The  actual  signal  as  it 
exists  in  space  is  complex  valued  and  has  a  single  sided  spectrum  for 
reasons  which  will  not  be  pursued  here. 

Fig  2  shows  the  modulus  of  the  complex  valued  spectrum  as  do  all  the 
subsequent  figures. 

Each  point  on  the  ground  reflects  a  replica  of  the  transmitted  chirp  and 
the  received  signal  at  any  instant  is  the  linear  sum  of  chirps  reflected  from 
points  illuminated  by  the  footprint  and  within  a  distance  corresponding  to  half 
a  pulse  width.  The  bandwidth  of  the  received  signal  is  therefore  the  same  as 
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that  of  a  single  chirp,  although,  of  course,  the  detail  of  the  spectrum  is 
greatly  different  and  is  noise-like  being  composed  of  chirps  with  random  phases 
and  displaced  in  time. 


The  actual  received  spectrum  is  sketched  in  Fig  7. 

The  receiver  coherently  mixes  the  signal  down  to  the  offset  video  frequency. 
This  offset  is  chosen  to  be  slightly  ^reater  than  AF/2  to  avoid  aliasing  the 
spectrum  and  is  one  quarter  of  t^e  subsequent  digitising  frequency  for  a  real 
valued  signal.  The  resulting  spectr  um  is  shew.-,  in  Fig  8. 

The  signal  is  still  t'-al  valued  at  this  point.  At  this  stage  the  signal 
could  be  operated  on  by  a.  quadrature  filter  to  give  the  analytic  signal  in  which 
case  the  negative  half  of  the  spectrum  disappears  and  then  f j  could  be  mixed 
down  to  zero  to  give  a  spectrum  symmetric  about  zero.  The  resulting  signal 
would  be  complex  valued.  For  the  moment,  however  a  real  valued  signal  will  still 
be  used. 

The  signal  is  next  digitised  at  a  sampling  frequency  f  (typically  exactly 
four  times  the  offset  video  fj  ).  The  process  of  sampling  the  signal  causes  the 
spectrum  to  be  periodically  .  epeated  at  a  frequency  of  f^  -  see  Fig  9.  Just  the 
envelope  of  the  spectrum  is  shown  in  Fig  9,  the  positive  and  negative  frequencies 
are  indicated  by  *  +  '  and 

When  range  compression  is  performed  in  the  frequency  domain  (as  it  nearly 
always  is)  using  a  discrete  Fourier  transform  (the  'fast*  Fourier  transform 
version)  then  the  spectrum  is  also  sampled.  If  N  samples  are  used  in  the 
Fourier  transform  the  spectrum  is  sampled  with  a  frequency  increment  of  f  /N 

5 

from  sample  to  sample,  see  Fig  10.  In  the  case  of  the  fast  Fourier  transform 
there  will  be  2n  samples,  typically  4096,  8192  or  16384  for  a  spacecraft  SAR. 
Notice  that  the  positive  frequencies  are  contained  in  samples  1  to  N/2  and  the 
negative  frequencies  in  samples  N/2  +  1  to  N  .  The  chirp  spectrum  is  similar 
and  thus  the  product  of  the  chirp  and  signal  spectra  is  similar. 

At  this  stage  the  negative  frequencies  are  eliminated  by  setting  samples 
N/2  +1  to  N  equal  to  zero,  and  on  inverse  transforming,  the  compressed 
analytic  signal  is  produced. 

6.2  The  azimuth  spectrum 

Previously  the  spectrum  of  just  one  pulse  was  considered.  Now,  the 
spectrum  of  a  large  number  of  pulses  is  discussed.  Fig  8  gives  the  spec t nun  of 
just  one  received  pulse  at  the  offset  video  frequency.  The  radar  periodically 
repeats  pulses  at  a  pulse  recurrence  frequency  of  v  Hz,  say.  The  effect  of  this 
is  to  turn  the  spectrum  in  Fig  8  into  a  line  spectrum.  The  distance  between  the 
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lines  is  v  Hz  .  The  resulting  spectrum  of  a  large  number  of  received  pulses  is 
shown  in  Fig  11,  This  should  not  be  confused  with  the  line  spectrum  already 
considered  in  Fig  10.  That  spectrum  arose  in  the  range  compression  process,  a 
process  independent  of  the  azimuth  direction. 

The  spectrum  shown  in  Fig  11  would  be  obtained  if  the  radar  platform  were 
stationary.  The  platform  moves,  of  course,  and  each  of  the  lines  in  the  spectrum 
therefore  rao\/es  due  to  the  Doppler  effect. 

Suppose  that  range  compression  is  performed  using  real  data.  The  result 
is  given  in  equation  (6''.  The  Doppler  ’phase’  term  is  a  cosine  function.  If 
range  compression  is  performed  with  an  analytic  signal  the  result  is  given  in 
equation  (7)  and  the  phase  term  is  an  exponential.  This  is  a  crucial  difference. 
Consider  just  one  line  in  the  spectrum  of  Fig  11.  If  the  real  ’phase’  term  is 
used  then  a  Doppler  shift  splits  the  line  into  two  lines  because  the  cosine  term 
implicitly  includes  both  positive  and  negative  frequencies  -  see  Fig  12.  In 
other  words  positive  and  negative  frequencies  are  indistinguishable,  ie  aliassed. 
For  the  analytic  signal  with  a  single  sided  spectrum  a  Doppler  shift  just  shifts 
each  line  as  shown  in  Fig  13.  There  is  only  one  line  and  posicive  and  negative 
frequencies  are  kept  separate. 

In  actual  fact  one  has  a  continuum  of  Doppler  frequencies  both  positive 
and  negative  for  a  broadside  looking  radar.  Each  line  is  then  broadened  into  a 
continuum  -  see  Fig  14. 

Note  that  the  radar  is  designed  in  such  a  way  that  the  pulse  recurrence 
frequency  is  high  enough  to  sample  the  Doppler  band  defined  by  the  antenna  beam- 
width  and  this  is  done  on  the  basis  of  an  analytic  signal. 

If  the  radar  is  a  pure  squint  mode  radar  with  only  positive  or  negative 
frequencies  then  the  aliassing  problem  would  not  exist.  One  could  then  either 
choose  a  high  enough  pulse  recurrence  frequency  or  a  very  carefully  chosen  prf 
so  that  the  Doppler  band  is  aliassed  into  a  wholly  positive  or  negative  band. 

This  would  undoubtedly  lead  to  problems  both  with  the  design  of  the  radar  and 
the  spacecraft  attitude  control. 

In  any  event  for  a  broadside  looking  radar  the  analytic  signal  is  always 
an  absolute  necessity  and  is  probably  also  necessary  for  practical  reasons  for 
a  pure  squint  mode  radar. 

6.3  Interpolation 

In  section  5.5  an  interpolation  technique  was  described  in  which  the  data 
to  be  interpolated  are  Fourier  transformed  via  a  DFT  and  then  inverse  transformed 
via  a  DFT  with  a  larger  number  of  samples.  In  using  this  technique  one  has  to  be 


37 


careful  to  treat  the  spectrum  correctly.  This  will  now  be  described,  with 
reference  to  interpolating  range  compressed  data. 

Consider  the  sampled  range  spectrum  shown  in  Fig  8.  It  corresponds  to  a 
Fourier  transformed  pulse.  Both  pulse  and  spectrum  are  sampled  and  periodically 
repeated.  The  period  of  the  spectrum  is  f  ,  the  sampling  frequency.  Suppose 
that  there  are  N  samples.  N/2  cover  the  negative  frequencies  and  N/2  the 
positive  frequencies,  see  Fig  9.  Increasing  the  number  of  samples  in  the  inverse 
transform  is  equivalent  to  sampling  at  a  higher  rate  and  the  spectrum  is  period¬ 
ically  repeated  at  a  longer  period  mf  ,  say,  where  m  is  an  integer  equal  to 
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some  power  of  2.  The  resulting  spectrum  is  shown  in  Fig  15.  In  order  to  achieve 
this  the  negative  frequencies  which  occupy  samples  N/2  +1  to  N  in  the 
original  spectrum  must  be  moved  up  to  samples  mN  -  N/2  +  1  and  mN  in  the  new 
spectrum  before  inverse  transforming  as  shown  in  Fig  15.  Everything  will  then 
be  correct  and  the  correct  interpolated  samples  will  be  produced. 

7  SIDELOBE  REDUCTION  -  APERTURE  AND  SPECTRUM  WEIGHTING 

Equation  (15)  gives  the  point  spread  function  in  both  range  and  azimuth  for 
the  ideal  matched  filtering  process.  The  modulus  of  the  point  spread  function  in 
both  directions  is  of  the  form  jsin  x/x|  .  This  function  is  plotted  in  Fig  16. 

It  will  be  observed  that  together  with  the  main  lobe  there  are  sidelobes.  These 
sidelobes  are  usually  considered  to  be  objectionable,  although  the  author  has 
processed  images  both  with  sidelobe  reduction  and  without  and  can  rarely  tell  the 
difference.  The  sidelobes  are  noticeable  in  the  case  of  very  bright  point-like 
scatterers. 

It  is  possible  to  reduce  the  level  of  these  sidelobes  but  only  at  the 
expense  of  broadening  the  main  lobe,  ie  reducing  the  resolution. 

If  the  correlation  process  (matched  filtering)  is  performed  in  the  time 
domain  then  an  integral  of  the  form: 


1  S<t)  | 


t0+T/2 
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exp  i2irxAF  —  dt 


t0-T/2 


(97) 


is  obtained.  If  the  correlation  process  is  performed  in  the  frequency  domain  an 
integral  of  a  similar  form  is  obtained  (this  is  the  inverse  Fourier  transform  of 
the  product  of  the  spectra  of  the  signal  and  chirp  replica  or  phase  history 
replica) . 
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where  the  time  centre  of  the  aperture  is  t^  ,  the  frequency  centre  is  fg  , 
the  time  width  is  T  ,  and  the  frequency  width  is  F  (this  is  equal  to  AF  but 
formula  (98)  has  been  written  as  it  has  in  order  to  exhibit  symmetry  with  (97). 

Sidelobe  reduction  is  obtained  by  weighting  the  integrals  in  (97)  and  (42). 
A  discussion  will  be  found  in  Ref  21,  section  3,4.2  p.780.  Such  weightit g  func¬ 
tions  are  usually  expressed  as  a  Fourier  series  with  a  period  of  either  T  or 
F  .  There  are  very  many  of  them.  The  ideal  function  (ideal  in  the  sense  that 
one  obtains  minimum  mainlobe  broadening  for  a  given  sidelobe  reduction)  is  the 
Dolph- Chebychef f  function  -  see  Ref  21,  section  3. 4. 2. 2  p.782.  In  Ref  16 
Papoulis  lists  a  number  of  different  functions  (section  7.3,  p.234). 

By  far  the  most  popular  weighting  function  is  a  'raised  cosine'.  This  will 
be  recognised  as  simply  the  first  two  terms  in  a  truncated  Fourier  series.  This 
it  also  known  as  a  Taylor  weighting  function  -  see  Ref  21. 

The  time  weighting  function  is: 


w(t)  =  1+28  cos 


pr(t-t0)-| 

T/2 


and  the  frequency  weighting  version  is: 


tQ  -T/2  <  t  <  tQ  +T/2  (99) 


w(f)  =  1+28  cos 
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fQ  -AF/2  <  f  <  fQ  +AF/2  (100) 


8  is  a  parameter  0  ^  8  <  £  which  specifies  the  sidelobe  levels, 
point  spread  function  is: 
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The  resulting 
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and  the  phase  part  of  the  psf  is  unaltered. 

2 

The  signal  energy  is  increased  by  a  factor  of  1+26  in  both  frequency 
and  time  domain  correlations  since 


I  “ 


2(t)dt  =  1  +  262  • 
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The  first  zero  in  the  point  spread  function  is  located  at  /l/(l  -2B)/AF  the 
position  of  the  other  zeros  are  unchanged.  The  relationship  between  3  and  the 
sidelobe  levels  must  be  obtained  numerically. 

An  alternative  to  the  Taylor  weighting  function  is  the  Chebyshev  weighting 
function,  as  pointed  out  above.  This  weighting  function  is  best  designed  by 
means  of  the  Remez  exchange  algorithm  -  see  Ref  23,  p.136.  The  sidelobe  levels 
are  specified  and  the  number  of  weights,  and  then  the  design  algorithm  gives  the 
weights.  Sidelobe  weighting  using  this  type  of  weighting  function  has  been  used 
in  Ref  9. 

8  RANGE  COMPRESSION 

The  various  aspects  of  the  range  compression  process  have  been  discussed  in 
previous  sections  and  now  in  this  section  they  are  collected  together. 

Range  compression  is  usually  performed  in  the  frequency  domain  using  the 
correlation  theorem  in  section  5.2  and  the  discrete  Fourier  transform  in  section 
5.3.  The  fast  Fourier  transform  (FFT)  version  of  the  disci«_t°  Fourier  transform 
(DFT)  is  invariably  employed  -  see  Ref  22  for  a  description  of  the  FFT.  The  FFT 
is  just  an  efficient  way  of  computing  a  DFT.  In  order  to  use  it  the  number  of 
samples  must  be  some  power  of  2,  ie  N  =  2n  . 

The  question  of  how  many  pulses  and  which  pulses  to  compress  is  considered 
in  the  sections  on  prefiltering  and  azimuth  compression. 

It  is  interesting  to  examine  the  difference  between  correlation  via  the 
time  domain  and  via  the  FFT.  The  number  of  multiplications  and  additions  for  the 
FFT  is  given  in  Ref  22.  If  there  are  N  complex  samples  in  the  signal  and  n 
complex  non-zero  samples  in  the  replica  then  correlation  via  the  FFT  requires 
approximately  N(l+log2N)  multiplications  and  2N  log2N  additions,  and  straight 
correlation  requires  nN  multiplications  and  nN  additions.  N  >  n  because  the 
replica  is  much  shorter  than  the  signal  for  obvious  reasons. 

As  an  example  consider  SEASAT.  If  one  chose  8192  real  range  samples  out  of 
a  total  of  13680,  ie  4096  complex  samples  for  the  analytic  signal,  and  1536  real 
value  samples  (768  complex)  for  the  chirp  replica  then  the  FFT  correlation 
requires  53248  multiplications  and  98304  additions  and  subtractions.  Straight 
correlation  requires  over  3  million  multiplications  and  3  million  additions. 
Correlation  via  the  FFT  is  therefore  much  faster. 

All  pulses  are  correlated  against  the  same  chirp  replica.  This  chirp 
replica  is  known  from  the  system  specifications.  Sometimes  it  has  to  be 
adjusted  slightly  in  order  to  achieve  the  best  possible  range  focussing. 


The  chirp  replica,  then,  is  sampled  at  the  same  rate  as  the  data.  The  chirp 
record  length  is  packed  out  with  zeros  to  make  it  the  same  length  as  each  pulse 
record  length  (=  2  ,  as  pointed  out  above).  The  packed  record  is  then  Fourier 
transformed  and  conjugated,  the  resulting  spectrum  is  then  weighted  by  the  side- 
lobe  weighting  function  given  in  section  7,  suitably  sampled.  Naturally  the 
spectral  width  of  the  weighting  function  AF  is  the  same  as  that  of  the  chirp 
spectrum.  The  weighted  chirp  spectrum  is  then  permanently  stored  and  used  to 
correlate  all  the  pulses.  Each  pulse  signal  is  taken  out  of  bulk  storage  and 
either  packed  with  zeros  to  2n  or  ♦  uncated  (this  is  more  usual).  Each 
pulse  is  then  Fourier  transformed,  multiplied  by  the  weighted  chirp  spectrum, 
and  inverse  Fourier  transformed.  This  process  is  continued  until  all  the  pulses 
required  to  process  the  image  have  been  compressed. 

The  reason  for  the  truncation  (or  packing)  of  the  pulse  samples  is  that  the 
total  number  of  samples  in  each  pulse  is  never,  in  practice,  an  exact  power  of  2. 
For  example  in  the  case  of  SEASAT  N  =  13680.  There  are  then  two  choices  and, 
usually,  considerations  of  total  storage  for  the  range  compressed  pulses  dictate 
that  each  pulse  be  truncated  to  8192  or  less. 

It  will  be  recalled  that  the  FFT  correlation  is  a  cyclic  correlation  - 
see  section  5.4,  and  this  will  now  be  examined.  The  signal  samples  can 
(schematically)  be  placed  at  equiangular  points  around  a  circle,  and  the  chirp 
replica  also,  see  Fig  17.  The  cyclic  convolution  then  implicitly  involves 
multiplying  each  signal  sample  with  the  corresponding  replica  sample,  adding, 
and  then  rotating  the  replica  relative  to  the  signal  by  one  sample  and  repeating 
the  process  N  times.  The  FFT  correlation  does  not  do  this  explicitly,  of 
course,  but  is  entirely  equivalent. 

If  the  signal  has  been  truncated  it  is  easily  seen  that  some  signal  has  been 
lost  for  a  time  of  T/2  at  each  truncation  point,  where  T  is  the  chirp  width. 

It  will  be  recalled,  that  due  to  the  finite  chirp  width,  the  signal  reflected  from 
each  point  scatter  does  not  extend  over  the  whole  signal  but  only  over  an  interval 
of  T  .  Hence  if  both  the  beginning  and  end  of  each  pulse  have  been  truncated  the 
signals  will  not  be  properly  correlated  over  cn  interval  of  ±T/2  around  the 
'join'  in  the  cyclic  correlation.  The  top  and  bottom  of  the  resulting  image  will 
then  be  spoilt  over  a  distance  corresponding  to  half  a  pulse  width.  It  is  there¬ 
fore  necessary  to  bear  this  in  mind  when  selecting  range  samples  to  cover  a  par¬ 
ticular  image  and  a  contingency  of  T/2  is  required  at  the  beginning  and  end  of 
each  pulse.  An  additional  contingency  is  required  to  cover  range  migration  and 
this  is  explained  further  in  section  13. 
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It  remains  to  be  shown  how  the  offset  video  frequency  f j  is  removed, 
see  equation  (7)  and  subsequent  remarks.  This  can  be  done  immediately  before 
inverse  Fourier  transforming  each  product  spectrum,  ie  after  each  signal  spectrum 
has  been  multiplied  by  the  replica  spectrum.  As  pointed  out  in  section  2.1  the 
offset  video  is  often  arranged  to  be  one  quarter  of  the  real  sampling  frequency. 
This  is  one  half  of  the  complex  sampling  frequency  (the  product  spectrum  at  this 
stage  is  that  of  the  analytic  signal).  If  the  offset  video  f^  is  exactly  half 
of  the  maximum  frequency  in  the  spectrum  then  to  translate  it  to  zero  the  spectrum 
is  simply  moved  by  one  half  period.  The  procedure,  then,  is  to  rotate  the 
product  spectrum  in  either  direction  by  N/2  before  inverse  Fourier  transforming. 

The  offset  video  can  also  be  removed  in  the  time  domain,  after  inverse 
transforming,  by  reversing  the  sign  of  alternate  samples.  Equation  (7)  gives  the 
phase  progression  due  to  the  offset  video  as  exp  i2irf,t  .  If  f„  is  the 
sampling  frequency  of  the  real  data,  fg/2  is  the  sampling  frequency  of  the 
complex  analytic  signal  and  the  time  increment  from  sample  to  sample  is  2/f^  . 
Hence  t  =  n(2/f  )  and  f,  =  f  /4  and  then  exp  i2uf1t  =  exp  iirn  =  (-l)n  . 
These  procedures  can  obviously  be  generalised  to  cope  with  the  situation  in  which 
the  offset  video  frequency  is  not  exactly  one  quarter  of  the  sampling  frequency. 

9  LOOK  FILTERING 

In  the  processing  of  an  image  the  Doppler  spectrum  is  usually  divided 
into  a  number  of  equal  parts  and  a  separate  image  is  processed  from  each  part. 

Thus  each  image  is  processed  from  a  sub-aperture.  The  radar  is  designed  to  have 
an  adequate  Doppler  bandwidth  to  permit  this  and  if  an  image  is  processed  using 
the  whole  Doppler  bandwidth  the  resulting  azimuth  resolution  is  many  times  finer 
than  the  range  resolution.  Each  image  processed  from  a  sub-aperture  is  detected 
{ie  the  modulus  is  formed,  or  an  intensity  image  formed)  and  then  the  images  are 
added  incoherently.  The  objective  is  to  produce  speckle  smoothing,  that  is,  an 
improvement  in  the  radiometric  resolution  at  the  expense  of  spatial  resolution. 

An  explanation  of  how  speckle  arises  will  be  found  in  Ref  24. 

Eacn  sub-aperture  is  generated  by  means  of  a  digital  band  pass  filter  -  the 
'look  filter'  -  centred  on  to  a  different  part  of  the  Doppler  spectrum.  The 
filtering  can  be  performed  before  the  azimuth  correlation  process  in  which  case 
it  is  known  as  prefiltering  or  alternatively,  during  the  azimuth  correlation 
process.  In  principle  it  could  also  be  carried  out  afterwards  but  this  is  never 
done  for  reasons  that  will  become  clear. 

Prefiltering  has  computational  advantages.  Suppose  that  the  aperture 
•.Doppler  spectrum)  is  divided  into  M  sub-apertures  each  of  width  AF  /M  where 
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AF  is  the  total  Doppler  speculum.  It  is  then  evident  that  since  we  have  an 
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analytic  signal  it  is  possible  to  subsample  the  Doppler  phase  and  reduce  the 
number  of  pulses  which  it  is  necessary  to  process.  Thus  for  M  sub-apertures 
it  is  only  necessary  to  process  every  Mth  pulse. 

9. 1  Design  of  the  prefilter 

The  first  requirement  of  the  prefilter  is  that  it  should  not  introduce 
spurious  phase  shifts  into  the  data,  since  this  will  have  serious  consequences 
for  the  azimuth  correlation  process.  A  linear  phase  response  (as  a  function  of 
frequency)  might  just  be  tolerable  since  this  would  just  shift  the  image  in  the 
azimuth  direction.  It  is,  however,  possible  to  design  a  digital  filter  which 
introduces  no  phase  shift  at  all.  This  is  one  of  the  advantages  of  operating 
with  digital  data;  a  zero  phase  shift  analogue  filter  is  impossible! 

There  are  basically  two  types  of  digital  filter,  the  finite  impulse  response 
(FIR)  and  the  infinite  impulse  response  (HR)  filters.  Only  the  first  type  can 
be  designed  to  have  zero  or  linear  phase  shifts.  It  is  possible  to  utilise  a 
filter  which  has  a  non  zero  phase  shift  by  passing  the  output  of  the  filter  back 
through  the  filter  in  reverse  time  sequence.  This  removes  the  phase  shifts  and 
is  called  two  pass  operation.  However,  this  filter  requires  all  the  output 
pulses  to  be  realised,  not  just  every  Mth  pulse.  In  addition  this  type  of  filter 
causes  difficulty  with  the  data  management.  There  are  similar  problems  with  the 
recursive  filter.  An  additional  problem  is  that  it  is  not  straightforward  to 
change  the  filter  characteristics  in  'mid  image'. 

By  far  the  most  straightforward  design  is  the  phase-shift  free  non  recursive 
filter  and  this  is  the  type  which  has  been  used  in  the  RAE  processor. 

FIR  filters  are  explained  in  detail  in  Ref  23. 

Consider  a  sequence  h(n)  -  ^  n  <  (“2~")  ^  °dd.  These  may  be 

regarded  as  the  filter  weights.  The  discrete  Fourier  transform  of  this  sequence 
is 

(N-I ) >2 

H(exp  jco)  =  h(n)  exp  -  icon  (103) 
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and  this  defines  a  finite  Fourier  series  period  2tt  and  gives  the  frequency 
response  of  a  non  recursive  filter.  If  there  is  no  phase  shift  then  obviously 
H(exp  ico)  is  real,  and  this  is  true  if  h(n)  =  h(-n)  and  then 
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This  argument  also  works  for  N  even  in  which  case 
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It  is  therefore  possible  to  design  a  filter  in  which  there  is  no  phase  shift. 

It  will  not  be  demonstrated  here,  but  the  case  N  even  leads  to  problems  because 
the  output  time  samples  occur  midway  between  the  input  samples.  Hence  filters 
used  for  prefiltering  have  an  odd  number  of  coefficients  and  have  symmetric 
coefficients. 

One  way  to  design  the  filter  is  the  'window  method'.  First  of  all  we 
decide  on  the  shape  of  the  required  response.  The  ideal  (from  one  point  of 
view)  is  the  rectangular  function  shown  in  Fig  18.  This  is  not  realistic, 
however,  because  it  requires  an  infinit  .lumber  of  Fourier  coefficients  and 
hence  an  infinite  sequence  of  input  samples.  The  resulting  Fourier  series  is 
truncated  therefore  and  then  the  coefficients  are: 
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where  v  =  pulse  recurrence  frequency,  u>c  =  cutoff  frequency  (rad/s) ,  and 
-(N-D/2  <n  <  (N-D/2  . 

Due  to  the  Gibb's  phenomenon  there  are  ripples  in  the  pass  band  and  in  the 
stop  band  and  the  amplitude  of  these  ripples  is  not  affected  by  the  order  of  the 
filter  N  .  These  ripples  may  be  smoothed  by  a  weighting  function  and  this 
function  works  in  exactly  the  same  way  as  the  sidelobe  weighting  function  in 
section  7.  Typically,  the  Hanning  weighting  function  is  used 


w(n)  =  1  +  cos  ~ 


(108) 
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Substituting  in  equation  (104)  then  gives  the  filter  response  function: 


H(exp  io))  = 


(1  +  cos  mr/N)  sin  nu)c/v  cos  nu/v 


(112) 


The  band  centre  can  be  shifted  by  weighting  by  an  exponential  weighting  function: 


h*(n)  =  h(n)  exp  - 


^irf^n 


(113) 


where  f^  is  the  shift. 

The  sub-sampling  and  weighting  procedure  is  then  described  by: 


(N+l)/2 

S'(m)  =  h(n)S(n  +  Mm) 

n=-(N-1 ) /2 


(114) 


M  is  the  data  reduction  factor  and  S(n),  n  =  1  ...  are  complex  range  compressed 
pulses.  S' (m)  are  the  prefiltered  pulses.  It  is  important  that  the  prefilter 
be  as  simple  as  possible  since  too  great  a  complexity  will  reduce  the  computa¬ 
tional  advantages  gained  from  data  reduction.  So  minimum  N  consistent  with 
sufficiently  low  ambiguities  is  chosen.  The  ambiguities  are  caused  by  non-zero 
response  in  the  stop  band,  ie  aliassing,  and  they  also  depend  on  the  data  reduc¬ 
tion  factor.  Response  curves  for  this  type  of  filter  are  given  in  Ref  11. 

Prefiltering  has  advantages  for  both  frequency  domain  azimuth  processing 
and  for  time  domain  processing.  The  reduction  of  the  data  rate  means  that  corner 
turning  (see  section  10)  can  be  done  on  smaller  arrays  with  consequent  time 
saving  and  the  forward  FFT  processes  a  smaller  amount  of  data.  The  savings  for 
time  domain  processing  can  be  considerable.  The  number  of  multiplications  in 
time  domain  processing  is  proportional  to  the  square  of  the  number  of  pulses, 
reducing  the  data  rate  by  4  (for  4  look  processing)  reduces  the  number  of  multi¬ 
plications  by  a  factor  of  16.  Likewise  the  number  of  multiplications  in  perform¬ 
ing  the  FFT  is  proportional  to  N  log2N  and  is  thus  reduced  by  a  factor  of  8. 

9.2  Frequency  domain  look  filtering 

Look  filtering  can  be  performed  in  the  frequency  domain,  and  this  approach 
is  commonly  followed  in  processors  which  perform  azimuth  correlation  in  the 
frequency  domain,  see  Ref  9  section  3.3.  An  azimuth  line  is  assembled  for 
azimuth  compression  (see  section  11  for  an  explanation  of  this  process)  and 
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Fourier  transformed  to  give  the  corresponding  spectrum.  This  spectrum  ther. 
covers  the  whole  Doppler  band.  If  M  looks  are  required  then  this  spectrum 
could  be  simply  split  up  into  M  non  zero  bands.  However,  the  opoortunity  can 
be  taken  here  to  introduce  sidelobe  weighting  at  the  same  time.  Simply  splitting 
up  the  spectrum  would  result  in  an  image  azimuth  point  spread  function  of 
'sin  x/x'  form  assuming  there  is  no  weighting  applied  to  either  the  Doppler 
spectrum  or  replica  (see  later  sections).  This  subject  is  further  explored  in 
the  sections  on  azimuth  compression. 

10  CORNER  TURNING 

After  range  compression  each  pulse  is  stored  on  a  direct  access  bulk 
storage  peripheral  -  invariably  a  disc  file  unit.  Such  disc  file  units  store  the 
data  serially  line  by  line.  Moreover  the  data  are  usually  transferred  to  and  from 
the  disc  file  in  the  form  of  blocks  of  typically  1024  16-bit  words.  These  blocks 
are  stored  at  random  on  the  disc  wherever  there  happens  to  be  space.  It  is  evi¬ 
dent  that  if  a  matrix  is  stored  row  by  row  and  it  is  required  to  access  the 
columns  of  the  matrix  then  very  many  blocks  will  have  to  be  retrieved  from  the 
disc  because  each  block  will  contain  few  column  elements.  Indeed  if  the  length 
of  each  row  in  the  matrix  exceeds  1024  words  (which  is  invariably  the  case  with 
SAR  data)  then  one  block  will  have  to  be  transferred  for  each  column  element. 

In  order  to  set  up  the  data  for  azimuth  compression  it  is  necessary  to  take 
samples  from  each  pulse,  i,e  to  read  along  the  stored  'columns'.  Unless  the  data 
are  organised  correctly  this  process  becomes  hopelessly  inefficient,  and  con¬ 
sequently  time  consuming. 

There  are  a  number  of  techniques  available  for  reorganising  the  data.  One 
of  the  simplest  and  more  efficient  methods  is  to  transpose  the  data  so  that  the 
stored  blocks  contain  the  data  in  serial  column  form  rather  than  serial  row  form. 
This  process  is  known  as  'corner  turning'.  This  transpose  operation  can  itselj. 
be  carried  out  in  a  number  of  ways.  The  technique  described  by  Eklundh^  is 
optimum  in  terms  of  the  number  of  operations  carried  out,  but  its  implementation 
depends  on  the  individual  computer  system  and  it  has  to  be  modified  to  handle  non 
square  matrices. 
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Corner  turning  is  not  strictly  necessary,  and  indeed  it  is  not  used  for 

example  in  the  RAE  processor.  It  is  used,  however,  in  both  the  Jet  Propulsion 

Laboratory  processor  in  the  USA^  and  the  processor  designed  by  Macdonald, 

9 

Dettwiler  and  Associates  in  Canada  .  One  alternative  strategem  is  to  write  special 
low  level  routines  to  transfer  data  on  and  off  disc  in  large  blocks,  taking  advan¬ 
tage  of  the  structure  of  the  disc  file  unit  to  incorporate  some  parallel  reading  and 
writing  to  and  from  the  columns  formed  by  the  tracks  on  the  stack  of  plattens 


\ 
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which  make  up  a  disc  file  unit.  The  size  of  the  blocks  is  limited  in  various  ways 
by  the  host  computer.  When  data  are  read  off  the  disc  for  azimuth  compression, 
then  this  is  done  in  large  (64  kbyte)  blocks.  The  shape  of  each  block,  ie  the 
row  and  column  lengths  is  chosen  to  maximise  the  number  of  new  image  points  which 
can  be  processed  (in  azimuth)  using  the  data  taking  into  account  the  range  migra¬ 
tion.  Then  all  the  image  points  which  can  be  processed  using  this  data,  and 
which  have  not  previously  been  processed,  are  processed. 

The  organisation  of  the  azimuth  data  for  both  frequency  and  time  domain 
processors  will  be  discussed  further  when  the  effects  of  range  migration  are 
examined . 

It  should  also  be  pointed  out  here  that  if  a  large  enough  random  access 
memory  is  available  to  store  all  the  range  compressed  data  to  process  a  given 
image  then  corner  turning  is  not  necessary.  At  the  time  of  writing  such  very 
large  memories  are  just  possible  but  they  are  expensive.  Matters  will,  however, 
change  and  in  a  few  years'  time  it  is  probable  that  corner  turning  will  become  a 
thing  of  the  past.  SAR  processor  designers  have  this  point  actively  in  mind  and 
future  designs  will  undoubtedly  feature  large  fast  random  access  memories. 

11  AZIMUTH  COMPRESSION  WITH  NO  RANGE  MIGRATION 

11.1  Frequency  domain  versus  time  domain 

As  pointed  out  in  section  2.2  azimuth  correlation  is  essentially  the  same 
process  as  range  correlation.  It  may  be  performed  either  in  the  time  domain  or 
in  the  frequency  domain  via  the  FFT.  For  a  small  number  of  samples  correlation 
via  the  time  domain  is  faster  than  via  the  FFT  due  to  the  inherent  speed 
advantage  of  the  FFT  method  being  outweighed  by  the  greater  time  penalty  of 
memory  read/ write  operations.  For  larger  N  ,  however,  the  FFT  method  is  much 
faster,  the  cross  over  point  depends  on  the  processor,  and  for  the  Floating  Point 
Systems  AP120B  for  example,  lies  between  64  and  128  points.  Satellite  SAR  azi¬ 
muth  processing  is  performed  on  data  sets  of  several  thousand  samples.  Hence, 
when  there  is  no  range  migration  there  are  clear  advantages  for  processing  in  the 
frequency  domain.  The  situation  when  there  is  range  migration  is  discussed  in 
section  14. 

11.2  Data  selection  and  packing 

Let  the  velocity  of  the  radar  platform  be  and  the  prf  v  .  Then  the 

distance  between  pulses  on  the  ground  is  V^/v  metres.  If  the  image  has  a  length 
L  in  the  azimuth  direction  then  Lv/L^  pulses  are  required,  ignoring,  for  the 
moment,  'end'  effects  due  to:  (a)  extra  pulses  at  each  end  of  the  record  being 
required  to  ensure  that  all  the  pulses  which  cover  the  image  points  at  the  ends 
of  the  azimuth  line  are  included,  and  (b)  to  cover  end  effects  in  the  correlation. 
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The  required  resolution  determines  the  Doppler  bandwidth  AF  -  see  equation  (15). 

“  2 

AF^  is  easily  shown  to  be  approximately  ( ie  to  first  order)  2V  T^/Aa^  where 
is  the  synthetic  aperture  time  interval,  a^  is  the  range  to  a  chosen  point 
at  t  =  0  ,  ie  broadside  on  in  this  case.  Hence  for  a  constant  AF^  ( ie  constant 
resolution)  must  increase  across  track  with  a^  .  This  happens  naturally 

because  the  projected  width  of  the  antenna  footprint  increases  linearly  with  a^ 
across  track,  see  Fig  19.  It  does  mean,  however,  that  in  order  to  maintain 
constant  azimuth  resolution  across  track  the  number  of  samples  in  the  azimuth 
data  {ie  the  number  of  probes)  must  be  increased  as  each  azimuth  line  is  pro¬ 
cessed.  In  addition  the  number  of  samples  in  the  azimuth  chirp  replica 
exp  -  i2irfQtp  (equation  (13))  must  also  be  increased.  Also  the  contingencies 
mentioned  above  to  cover  the  ends  of  the  azimuth  line,  are  also  range  dependent. 


Consider,  for  the  moment,  the  processing  of  the  whole  azimuth  spectrum  to 
obtain  a  maximum  resolution  one  look  image.  If  the  azimuth  length  of  the 
required  image  is  L  ,  the  prf  v  and  the  radar  velocity  V  then  evidently 
the  number  of  pulses  required  is  Lv/V  .  The  number  of  samples  in  the  Doppler 
replica  against  which  the  signal  is  correlated  is  given  by  the  aperture  time 
T^  =  AF^AaQ/2Vp  and  the  prf,  v 


That  is  the  number  of  samples  =  vT^ 


The  contingencies  should  ideally  be  added  to  the  longer  of  the  two  records, 

ie  if  the  number  of  pulses  n  >  n^  ,  the  number  of  replica  samples  then  the 

contingency  (=  n^/2  for  each  of  the  two  contingencies  (1)  and  (2)  given  above, 

ie  n  ) ,  should  be  added  to  each  end  of  the  n  record. 
r  P 

In  order  to  perform  the  correlation  via  an  FFT  the  number  of  samples 

should  be  2°  and  one  has  the  option  of  either  packing  the  record  length 

with  zeros  or  choosing  a  record  length  (including  contingency)  equal  to  2°  . 

In  fact  both  the  signal  record  (azimuth  line)  and  azimuth  replica  expand  as 

each  azimuth  line  is  processed  moving  across  track  in  order  to  maintain 

resolution. 


If  one  is  prepared  to  accept  a  variation  in  azimuth  resolution  across 
track  then  constant  length  records  are  acceptable.  In  this  case  it  is  only 
necessary  to  pack  one  record  (probably  the  replica)  whereas  in  the  former  case 
it  is  necessary  to  pack  both.  However  variable  length  records  do  not  constitute 
a  major  problem. 

11.3  The  replica 

The  replica  is  exp  -  ^rTf^t^  with  t^  given  by  equation  (13)  except  that 
in  this  section  it  is  assumed  that  a^  =  0  ,  ie  there  is  no  range  walk.  The 
replica  is  then  the  same  for  all  points  along  a  given  azimuth  line  at  constant 
range  sample  number,  the  only  difference  being  a  shift  in  the  time  axis.  This 
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translates  into  a  phase  shift  in  the  spectrum  and  therefore  all  points  con¬ 
tribute  the  same  spectrum  (scaled  by  the  scattering  cross  sections)  with  differ¬ 
ent  phase  shifts.  Thus  all  points  may  be  correlated  by  the  same  replica 
spectrum  and,  on  inverse  transforming,  these  phase  shifts  translate  back  into 
time  shifts  and  each  focussed  point  is  returned  to  its  correct  sequence  in 
the  image  line. 

The  range  history  is  computed  for  each  azimuth  line,  but  is  constant  for 
that  line  once  computed.  It  is  emphasised  that  in  this  section  it  is  assumed 
that  there  is  no  range  migration  effect  (which  is  valid  for  most  aircraft  SARs 
for  instance).  The  replica  is  computed  at  time  increments  equal  to  the  time 
between  each  pulse  and  the  number  of  samples  corresponds  to  the  number  of  pulses 
in  the  aperture,  or  sub  aperture. 

11.4  The  correlation  process 

This  is  identical  to  the  range  compression  process  described  in  section  8. 
The  Doppler  signal  (the  azimuth  line)  is  Fourier  transformed  via  an  FFT.  The 
range  history  is  computed  and  Fourier  transformed,  complex  conjugated  and  the  two 
spectra  are  multiplied  together,  multiplied  by  the  sidelobe  reduction  weights 
(section  7),  inverse  transformed,  and  the  image  line  stored.  The  end  image 
points  in  the  record  corresponding  to  half  the  replica  length  are  rejected  since 
they  are  not  properly  correlated. 

11.5  Multilook  processing 

So  far  'one  look'  processing  with  the  entire  Doppler  spectrum  has  been 
assumed.  In  practice  multilook  images  are  usually  required  and  each  separate 
look  is  produced  by  processing  a  sub-aperture.  This  can  be  done  by  processing 
prefiltered  data  (see  section  9)  or  by  selecting  a  sub-aperture  during  the 
correlation  process  (as  is  done  in  Ref  9  for  instance) . 

When  processing  prefiltered  data  the  processing  proceeds  exactly  as  before 

except  that  the  effective  prf  is  now  reduced  by  the  data  reduction  factor  (which 

equals  the  number  of  looks).  Thus  to  process  a  given  number  of  image  points 

requires  fewer  pulses.  Naturally,  each  look  is  processed  from  a  different  set 

of  prefiltered  data  so  that  the  total  number  of  pulses  processed  remains  the 

same.  However  overall  processing  time  is  shorter  because  it  depends  upon  the 

number  of  pulses  in  a  non-linear  way.  The  effective  prf  is  reduced  by  the  data 

reduction  factor  M  and  there  are  fewer  pulses  in  the  aperture,  but  since  the 

(sub)  aperture  is  now  reduced  in  width  also  by  a  factor  M  the  total  number  of 

2 

pulses  in  the  aperture  is  reduced  by  M  (16  for  a  4  look  system).  Hence  if 
there  are  approximately  4000  pulses  in  the  full  aperture  (for  SEASAT)  there  are 
only  250  effective  pulses  in  one  quarter  aperture  (approximately).  For  an  eight 
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look  system  there  would  only  be  63  or  so  and  time  domain  correlation  would  begin 
to  become  a  possibility! 

Hence  in  the  construction  of  the  replica  the  range  history  is  computed  at 
intervals  of  At/M  where  At  =  1/v  and  only  for  those  samples  within  an  interval 
of  T/M  where  T  =  full  aperture  width.  Note  that  each  sub-replica  is  computed 
using  exactly  the  same  parameters,  the  sub-replicas  only  differ  in  the  range  of 
time  over  which  they  are  computed.  This  automatically  ensures  that  the  azimuth 
image  lines  computed  from  each  sub-aperture  are  correctly  registered.  The  process¬ 
ing  then  proceeds  as  before.  The  sub-replica  is  packed  with  zeros  to  make  ic 
the  same  length  as  the  azimuth  line  (in  terms  of  the  numbers  of  samples),  Fourier 
transformed,  multiplied  by  the  sidelobe  reduction  weights,  multiplied  by  the  com¬ 
plex  conjugate  of  the  Fourier  transform  of  the  azimuth  line  and  then  inverse 
transformed.  This  is  repeated  for  the  other  sub-apertures.  The  imperfectly 
correlated  image  points  »*■  the  beginning  and  end  of  the  image  lines  are  rejected. 
The  modulus  of  each  of  ^ne  complex  valued  azimuth  lines  is  taken  and  they  are 
added  incoherently  to  produce  a  smoothed  image  line. 

An  alternative  strategy  to  prefiltering  the  data  is  to  generate  the  looks 
during  the  correlation  process.  Again,  the  processing  proceeds  much  as  in  the 
full  aperture  case.  The  pulses  covering  the  required  azimuth  line  (all  of  them) 
are  selected  and  a  contingency  added  at  each  end.  They  are  packed  to  2n 
and  Fourier  transformed.  The  spectrum  is  then  split  up  into  a  number  of 
separate  contiguous  sections.  These  are  the  looks.  The  replicas  are  constructed 
exactly  as  in  the  case  of  the  prefiltered  data  and  the  correlation  is  then 
performed  in  the  same  way.  This  technique  is  used  in  Ref  9:  it  will  be  noted 
that  the  sidelobe  reduction  filter  used  in  Ref  9  was  of  a  different  kind  from  that 
suggested  here,  having  a  linear  phase  equal  ripple  Chebyshev  response,  see 
section  9.2. 


11.6  Azimuth  correlation  point  by  point 

This  technique  is  very  slow  in  comparison  with  the  block  processing 
inherent  in  the  FFT  correlation.  It  is  not  limited  to  a  fixed  pixel  spacing  of 

v/afa  . 

P  A 

The  pulses  covering  the  target  point  are  selected  and  the  azimuth  phase 
replica  computed  for  the  image  point  Tor  which  the  data  are  being  correlated. 

This  is  a  one  point  correlation  and  the  time  for  which  the  replica  is  computed 
corresponds  to  the  time  coordinate  of  the  imaged  point  in  the  image.  The  replica 
samples  are  complex  conjugated  and  sidelobe  weighted  and  the  resulted  weighted 
samples  are  each  multiplied  by  the  corresponding  azimuth  signal  samples  (of 
which  there  art  an  equal  number) .  The  resulting  samples  are  then  summed  and 
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this  corresponds  to  the  correlated  image  point.  The  processor  then  moves  on  to 
the  next  point  (which  is  arbitrary)  and  repeats  the  process. 

12  MAPPING  AND  INTERPOLATION 

The  previously  explained  range  and  azimuth  correlation  techniques  produce 
an  image  in  slant  range  at  fixed  pixel  spacings  of  c/2AFR  and  in  azimuth  at  a 
fixed  pixel  spacing  of  V^/AF^  where  c  and  are  the  velocities  of  light 

and  the  radat  platform,  and  AFR,  Af^  are  the  video  (chirp)  bandwidth  and 
Doppler  (sub-aperture)  bandwidth.  See  equation  (15)  and  subsequent  comments. 

There  are  two  difficulties  here.  First,  image  users  invariably  require 
the  image  to  be  an  orthographic  map  projection  and  not  mapped  into  slant  range. 
Second,  arbitrary  pixel  spacing  is  desirable  not  the  least  because  the  image 
detection  process  (ie  finding  the  modulus)  increases  the  effective  bandwidth 
when  the  data  is  converted  from  complex  form  to  real  form.  It  is  therefore 
necessary,  in  principle,  to  increase  the  image  sampling  rate. 

The  explicit  way  to  do  this  is  simply  to  interpolate  between  image  sample 
points  using,  for  example,  cubic  splines  (see  Ref  26  p.349)  or  by  convolving 
with  an  interpolating  function  which,  in  the  case  of  a  SAR  image,  would  be  the 
point  spread  function.  In  fact  it  is  far  better  to  interpolate  an  image  before 
detection,  ie  by  interpolating  the  complex  valued  image.  At  least,  then,  the 
image  is  --rrectly  sampled  before  the  interpolation  begins!  There  are  several 
possibilities  and  one  of  the  most  practical  is  to  weight  the  point  spread 
function  of  the  system  by  several  samples  on  either  side  of  the  point  at  which 
the  image  is  required  to  be  interpolated,  add  the  point  spread  functions  and 
evaluate.  This  is  just  a  limited  convolution  over  a  few  samples. 

One  technique  which  is  in  use  at  RAE  in  the  range  direction  is  that 
explained  in  section  5.5.  After  range  compression  an  inverse  FFT  is  performed 
of  k  times  the  initial  length.  This  automatically  generates  k  -  1  inter¬ 
mediate  samples  between  the  original  samples.  The  subsequent  data  selection 
then  chooses  range  samples  nearest  to  the  desired  ones.  This  technique  is  a 
very  fast  way  of  performing  interpolation  but  it  requires  a  lot  more  memory.  It 
is  a  powerful  technique  when  used  in  conjunction  with  data  which  has  range  walk 
because  these  intermediate  samples  reduce  ghosting  (see  section  15)  and  in 
addition  the  range  walk  implicitly  interpolates  the  data  (see  section  13), 

Most  interpolation  methods  give  trouble  at  the  ends  of  a  data  set. 

In  principle  given  an  infinite  data  set  sampled  at  equal  intervals  the 
original  bandlimited  fraction  can  be  completely  reconstructed  from  its 
samples  (see  Ref  16,  p  141),  simply  by  convolving  the  data  with  a  'sin  x/x' 


function.  In  practice  one  has  a  finite  data  set.  It  can  be  shorn  that 
convolving  such  a  finite  data  set  with  a  'sin  x/x'  function  is  equivalent  to 
fitting  a  polynomial  of  degree  n  -  1  to  the  n  samples.  (Actually  this  is 
Lagrangian  interpolation,  see  Ref  26  p.235).  The  difficulty  is  that  such  poly¬ 
nomials  'wiggle'  and  although  the  interpolation  may  be  smooth  and  accurate  in  the 
middle  of  the  data  set,  at  the  edges  it  could  be  a  long  way  out.  Fortunately  the 
ends  of  the  data  set  are  not  properly  correlated  anyway  due  to  end  effects  on  the 
circular  convolution  and  in  practice,  bearing  in  mind  that  the  end  points  of  the 
correlated  data  are  discarded  anyway,  this  method  seems  to  work  very  well.  It 
will  be  recognised  that  inverse  Fourier  transforming  with  a  larger  interval  than 
the  forward  transform  effectively  convolves  with  a  'sin  x/x'  function  and  evaluates 
in  between  the  original  samples.  The  situation  is  also  helped  by  sidelobe 
weighting  applied  to  range  compression  which  keeps  down  the  sidelobes  in  the 
point  spread  function  and  hence  the  wiggles  in  the  interpolating  polynomial. 

In  the  azimuth  direction  one  can  also  interpolate  implicitly  by  performing 
the-azimuth  correlation  point  by  point. 

1 3  RANGE  MIGRATION 

The  Doppler  effect  is  a  result  of  the  changing  phase  between  the  radar  and 
the  target  as  the  radar  moves  past.  The  phase  <f>(t)  is  given  by 

<j>(t)  =  (a0  +  ajt  +  a2t2  +  ...)  (115) 

a^,  aj  and  a2  are  defined  after  equation  (10). 

The  changing  phase  is  a  result  of  the  change  in  slant  range  between 
radar  and  target.  If  the  change  in  range  is  smaller  than  approximately  one  half 
of  the  range  resolution  length  over  the  whole  synthetic  aperture  then  the  samples 
used  in  the  azimuth  compression  process  (the  azimuth  line)  can  be  taker,  from  the 
same  sample  number  (or  range  gate)  from  each  pulse.  The  factor  of  one  half  is, 
introduced  by  the  double  path  length. 

On  the  other  hand  if  the  change  in  slant  range  is  greater  than  one  half 
resolution  cell  then  the  azimuth  samples  used  in  azimuth  compression  have  to  be 
taken  from  different  range  gates  in  each  pulse.  This  effect  is  known  as  range 
migration  and  has  already  been  discussed  from  a  different  point  of  view  in 
section  2,3.  Evidently  the  lccus  of  a  point  through  the  range  gates  (sample 
space)  is  very  nearly  parabolic  since  the  Doppler  shift  is  very  nearly  a  linear 
function  of  time  with  a  constant  offset.  The  shape  of  the  parabola  is  given  by 
equation  (115)  and  contains  a  linear  component  (the  aj  term)  and  a  quadratic 
component  (the  a2  term).  In  this  Report  the  linear  component  is  called  range 


walk  and  the  quadratic  component  (the  a*,  term)  is  called  range  curvature. 
Together  they  constitute  range  migration. 

The  behaviour  of  the  a^  and  a9  components  were  investigated  in  sections 


4.3  and  4.4.  Equation  (66)  shows  that  the  range  walk  component  is  caused  by 
orbit  eccentricity,  squint  and  Earth  rotation.  The  rotation  term  is  dependent  on 
the  latitude  because  the  Earth's  surface  velocity  is  greatest  at  the  equator  and 
smallest  at  the  'top'  of  the  orbit  (nearest  the  poles)  because  of  the  much  smaller 
projected  Earth  radius.  In  addition  the  included  angle  between  the  swath  velocity 
vector  and  the  Earth's  surface  velocity  varies,  being  equal  to  the  orbit  inclina¬ 
tion  angle  at  the  equator  (if  the  squint  is  zero)  and  zero  nearest  the  poles. 

The  range  walk  is  hence  a  function  of  latitude  and  is  not  constant  in  the  along 


track  direction,  although  it  can  be  considered  so  for  a  small  enough  interval. 
The  range  walk  is  also  a  function  of  a^  and  of  0q  ,  ie  across  track,  since 
the  grazing  angle  varies  and  hence  the  resolved  component  of  velocity. 

The  quadratic  component  is  also  a  function  of  latitude  although  a  much 


weaker  function  than  the  linear  component. 

Hence  in  the  case  of  an  orbital  SAR  the  range  migration  consists  of  two 
components  (1)  a  linear  component  (range  walk),  and  (2)  a  quadratic  component 
(range  curvature).  The  first  is  a  result  of  the  rotating  Earth,  non  circular 
orbit  and  squint  and  the  second  is  primarily  a  result  of  the  changing  geometry. 
Note  that  curvature  could  also  be  caused  by  higher  order  terms  a^,  a^,  etc. 
These  have  been  ignored  here;  the  cubic  curvature  term  a^  was  investigated  in 
section  4.4  and  for  the  purposes  of  range  migration  correction  it  is  negligible. 

The  situation  is  made  clearer  by  Fig  20  (which  has  been  adapted  from 
Ref  8  Fig  1).  In  Fig  8  Aw  is  the  total  range  walk  and  Ac  is  the  range 
curvature. 

The  range  to  a  target,  is  given  by: 

a(t)  =  aQ  +  ajt  +  a2t  .  (116) 

For  the  example  discussed  in  sections  4.3  and  4.4  a.  =  65.66  m/s  and 
a2  =  28.09  m/s  .  For  SEASAT  the  distance  between  range  gates  (samples)  is 
about  7.89  m.  Hence  over  an  interval  of  2  s,  corresponding  to  the  maximum 
synthetic  aperture,  the  total  linear  range  walk  Aw  is  about  65.66  x  2  *  2/7.89 
range  samples  =  33.3.  The  extra  factor  of  2  is,  again,  due  to  the  double  path 
length.  In  the  case  of  the  curvature  we  take  half  of  the  interval  or  1  s  and 
this  gives  Ac  =  7.1  range  gates.  Note  that  the  range  walk  can  be  much  greater 
nearer  the  equator. 
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The  range  migration  problem  is  the  dominant  problem  in  the  design  of  digi¬ 
tal  SAR  processors  for  satellite  data;  and  if  it  can  be  ignored  it  very  greatly 
simplifies  the  processor  design.  For  the  future  it  may  be  that  new  technology 
will  permit  different  satellite  systems  in  which  the  effect  can  be  ignored,  if  so 
it  will  be  possible  to  much  simplify  the  processor  design. 

14  AZIMUTH  CORRELATION  WITH  RANGE  MIGRATION 

Azimuth  correlation  with  range  migration  can  either  be  performed  in  the 
frequency  domain  or  in  the  time  domain.  The  frequency  domain  approach  is 
examined  first  as  this  is  the  most  usual  technique.  The  method  described  in  this 
section  is  that  used  in  the  JPL  processor,  see  Refs  7  and  8,  it  is  also  employed 
in  the  MDA  SEASAT  processor,  see  Ref  9. 

14.1  Azimuth  correlation  in  the  frequency  domain 

The  range  migration  curve  shown  in  Fig  20  is  shaped  like  a  banana  and  is 
therefore  commonly  known  at  RAE  as  the  data  'banana'.  That  term,  being  short  and 
succinct,  is  used  here.  The  box  of  data  containing  the  banana,  that  is  the  rec¬ 
tangular  box  of  data  within  which  the  banana  can  be  fitted  is  called  the  'banana 
box' . 

The  first  essential  idea  behind  the  simultaneous  processing  of  a  number  of 
azimuth  image  points  via  the  FFT  is  that  the  banana  only  slowly  varies  along  the 
azimuth  direction.  That  is  a  number  of  consecutive  targets  all  having  the  same 
a^  have  nearly  the  same  banana.  Obviously  the  bananas  are  shifted  in  azimuth 
relative  to  each  other  but  the  shape  of  each  banana  is  approximately  the  same. 
This  is  shoT-?n  in  Fig  21. 

The  second  essential  idea  is  that  there  is  a  one  to  one  correspondence 
between  the  azimuth  coordinate  (in  time)  and  the  corresponding  Doppler  frequency 
of  a  target.  This  is  obvious  from  equations  (8)  and  (9).  Therefore  if  a  number 
of  bananas  are  Fourier  transformed,  the  same  points  in  each  banana  map  onto  the 
same  point  in  Doppler  frequency/range  space.  This  is  shown  in  Fig  22.  One  can 
then  define  a  spectrum  banana  and  a  corresponding  spectrum  banana  box. 

Hence  if  one  has  a  box  of  data  containing  many  bananas  corresponding  to 
many  targets  all  with  the  same  a^  ,  and  these  bananas  are  Fourier  transformed 
the  spectra  map  on  top  of  each  other  with  relative  phases  corresponding  to  their 
relative  displacements.  The  spectrum  banana  is  then  formed.  Evidently  other 
data  bananas  corresponding  to  different  values  of  ap  are  Fourier  transformed  in 
the  same  way  and  mapped  one  to  one  into  the  spectrum  box. 

In  practice  one  does  not  wish  to  Fourier  transform  along  a  banana  but 
across  an  azimuth  line  at  constant  range.  If  this  is  done  it  will  be  seen  from 
the  above  arguments  that  the  spectrum  bananas  are  automatically  created  in  the 
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spectrum/ range  space  due  to  the  one  to  one  correspondence  between  offset  along 
the  banana,  range  and  Doppler  frequency. 

It  is  straightforward  to  identify  a  spectrum  banana  corresponding  to  a 


group  of  data  bananas  because  the  spectrum  banana  is  simply  a  data  banana 


stretched  out  in  the  azimuth  direction  to  fit  the  spectrum  width. 


Tnis  elegant  technique  will  only  work  if  all  the  data  bananas  have  the  same 
shape  and  herein  lies  a  difficulty.  This  requirement  sets  the  maximum  azimuth 
distance  which  can  be  correlated  simultaneously.  The  data  are  thus  divided  up 
into  sections.  Each  section  has  a  set  of  different  bananas  with  different  range 
walk  gradient.  The  sections  must  overlap  in  order  to  ensure  that  all  target 
points  are  properly  correlated  (end  effects). 

Once  the  box  of  data  has  been  selected  and  Fourier  transformed  line  by  line 
the  spectrum  bananas  are  obtained  by  piecewise  approximation.  A  section  of 
spect-rum  is  taken  from  the  appropriate  part  of  each  line.  Each  section  of 
spectrum  is  then  coherently  assembled  to  form  a  complete  spectrum  banana  -  see 
Fig  21.  The  appropriate  section  of  spectrum  free  each  line  is  selected  by  taking 
that  part  of  each  line  which  has  within  one  half  range  gate  of  the  required 
spectrum  banana.  Note  that  all  of  this  is  possible  because  of  the  one  to  one 
correspondences  noted  above  and  the  linear  nature  of  the  Fourier  transform. 

After  the  spectrum  has  been  assembled  it  is  operated  on  in  exactly  the  same 
way  as  in  section  11,  to  which  the  reader  is  now  referred.  Briefly,  the  'looks' 
are  generated,  unless  the  data  is  prefiltered,  and  the  replica  spectrum  is  generated 
via  the  range  history  and  Fourier  transformed,  conjugated  he  spectra  multiplied 
and  inverse  transformed.  The  resulting  image  line  is,  of  ;urse,  mapped  on  to 
the  image  along  a  straight  line  at  tie  corresponding  value  of  . 

14.2  Pre-skewing 

It  is  evident  that  the  number  of  spectrum  segments  required  is  largely 
dependent  on  the  range  walk.  An  aperture  may  have  a  walk  of  several  tens  of 
range  gates  and  there  will  then  be  several  tens  oi  su'-ccrum  segments.  It  is 
desirable  to  decrease  the  number  of  segments  to  improve  the  epeed  of  the 
processing.  One  way  to  do  this  is  to  pre-skew  the  data  by  rearranging  the  rang., 
gates  relative  to  each  other  -  see  Fig  24.  The  rate  of  displacement  of  the  range 
gates  is  chosen  to  be  equal  to  the  ratine  walk.  It  is  then  only  necessary  to 
produce  a  segmented  approximation  to  the  range  curvature  instead  of  both  the  walk 
and  the  cur\.  ture.  This  greatly  reduces  the  number  of  segments. 

This  technique,  however,  has  a  severe  disadvantage.  The  range  walk  is 
latitude  dependent  as  pointed  out  in  section  13  and  the  data  must  therefore  be 
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divided  into  sections  as  in  section  14.1.  Unfortunately  there  are  difficulties 
in  fitting  the  sections  together  with  the  pre-skewing  method  because  of  the 
different  range  walk  gradient  from  section  to  section.  There  are  small  gaps 
between  sections  or  alternatively  non  registering  overlaps  and  extra  manipula¬ 
tion  is  required  to  make  the  sections  fit  together. 

This  problem  does  not  arise  if  the  range  migration  correction  is  done 
entirely  in  the  frequency  domain  -  see  section  14.1. 


» 
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14.3  Azimuth  correlation  in  the  time  domain 

This  is  the  technique  used  in  the  RAE  processor^.  It  is  the  most  straight¬ 
forward  way  to  correlate  data  which  has  severe  range  migration.  Each  image  point 
is  correlated  separately.  The  azimuth  data  is  selected  for  the  image  point  of 
interest  by  computing  the  data  banana  and  then  selecting  each  range  sample  lying 
within  a  specified  fraction  of  a  range  gate  of  the  required  banana.  If  the 
range  compressed  data  is  not  interpolated  this  fraction  is  one  half,  see  Fig  25. 
Thus  in  contrast  to  the  frequency  domain  correlation,  it  is  the  data  banana 
which  is  piecewise  approximated.  When  the  data  banana  has  been  assembled  the 
range  history  of  the  point  is  computed  at  the  time  corresponding  to  the  relevant 
time  coordinate  in  the  image.  The  banana  samples  and  the  complex  conjugate  of 
the  replica  samples  are  then  multiplied  together  and  this  gives  the  single  point 
correlation  for  this  image  point.  The  process  is  repeated  for  the  next  image 
point  with  a  new  banana  and  a  new  replica.  The  next  image  point  can  be  spaced 
completely  arbitrarily  from  the  preceding  point  both  in  azimuth  and  range.  The 
range  migration  effect  helps  to  interpolate  the  data  implicitly. 

It  will  be  seen  that  this  technique  is  very  precise  and  involves  no 
fundamental  approximation.  There  is  the  piecewise  approximation  of  the  data 
banana  but  the  effects  of  this  can  be  mitigated  by  interpolating  the  range 
samples  (using  the  Fourier  transform  interpolating  technique  described  in 
section  5.5) . 

Time  domain  correlation  is  slow.  However,  there  are  considerable  advan¬ 
tages.  First,  the  sample  spacing  of  the  image  is  completely  arbitrary.  One 
result  of  this  is  that  the  image  can  be  mapped  in  any  way  without  interpolation. 
Second,  the  image  can  be  rotated,  the  x  and  y  axes. of  the  image  do  not 
necessarily  have  to  align  with  the  range  and  azimuth  directions.  Third,  there 
are  no  awkward  problems  involved  in  fitting  sections  together.  Above  all 
tbo  image  is  precisely  correlated.  All  these  advantages  together  produce  images 
of  very  high  quality. 
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15  GHOST  IMAGES 

uhosts  are  faint  replicas  of  the  main  image  displaced  in  the  azimuth 
direction  and  appear  on  either  side  of  the  main  image.  Fig  26  shows  the  image 
of  a  point  target  (it  is  the  same  image  as  Fig  1 ,  but  without  ghost  suppression) 
and  its  attendant  ghosts.  There  are  usually  many  ghosts  visible  and  they  are 
spaced  out  from  the  main  image  at  regular  intervals.  The  position  of  the  ghosts 
depends  upon  the  sub-aperture  position  and  Fig  26  shows  three  images  from  three 
adjacent  sub-apertures  added.  The  ghosts  then  appear  in  sets  of  three.  Fig  26 
only  shows  clearly  the  fundamental  ghosts,  but  in  fact  there  are  other  groups  of 
three  spaced  out  further  from  the  central  point  which  are  not  easily  visible. 
These  ghost  artifacts  are  produced  during  the  azimuth  correlation  process  and 
are  caused  by  the  piecewise  approximation  of  the  spectrum  banana  in  frequency 
domain  processing  or  the  piecewise  approximation  of  the  data  banana  in  the  time 
domain,  see  section  14  for  an  explanation  of  these  approximations.  Ghosts  are 
produced  by  both  frequency  domain  and  time  domain  processors. 

The  artifact  is  briefly  analysed  in  the  section  by  considering  the  image 
of  a  single  point  on  an  absorbing  background.  Fig  27  shows  the  effect  of  the 
piecewise  approximation  in  the  time  domain,  and  Fig  28  shows  the  corresponding 
effect  on  the  frequency  domain.  In  both  cases  the  data  cycles  over  the  top  of 
the  range  point  spread  function  and  hence  the  data  assembled  for  azimuth  com¬ 
pression  is  amplitude  modulated.  There  is  also  a  phase  modulation,  but  this  will 
be  ignored  here.  The  peak  to  peak  ripple  amplitude  for  a  'sin  x/x’  point  spread 
function  in  the  range  direction  and  range  sample  spacing  equal  to  the  distance 
between  the  zeroes  of  the  point  spread  function  (=c/2AF  )  is  0.363  of  the  ampli- 
tude  of  the  psf,  ie  -8.8  dB.  The  fundamental  frequency  of  the  modulation  is 
given  by  the  rate  at  which  the  data  cycles  over  the  point  spread  function.  The 
analysis  which  follows  is  applicable  to  a  time  domain  processor,  it  can,  however, 
be  easily  modified  to  make  it  apply  to  the  frequency  domain.  The  range  is  given 
by 

a(t)  =  aQ  +  ajt  +  a2t2  .  (1  17) 


The  distance  between  range  samples  is  c/2AF  and  so  the  frequency  is  given  by 
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The  amplitude  modulation  produces  coherent  sidebands  in  the  time  domain  data. 


The  Doppler  information  is  contained  in  the  sidebands  and  these  are  also  corre¬ 
lated  but  being  shifted  in  frequency  from  the  main  Doppler  band  they  are 
shifted  in  the  azimuth  direction  in  the  image.  These  coherent  sidebands  are, 
then,  the  source  of  the  ghosts  in  the  image.  Unfortunately  in  order  to  give 
a  satisfactory  analysis  of  the  ghosts  it  is  necessary  to  include  second-order 
effects  to  explain  the  observed  shift  of  the  ghost  positions  with  sub-aperture 
position.  Such  an  analysis  can  be  performed  but  is  too  long  to  present  here 
(see  Ref  10  for  additional  information). 

The  most  elegant  way  to  reduce  the  effect  of  ghosting  is  to  interpolate 
the  range  compressed  data.  This  reduces  the  period  of  cycling  on  the  piecewise 
approximation,  te  increases  the  frequency.  As  a  result  the  cycling  occurs 
over  a  smaller  segment  of  the  point  spread  function  and  this  greatly  reduces  the 
amplitude  of  the  ripple.,  The  type  of  interpolation  referred  to  here  is  simply 
the  generation  of  intermediate  samples  between  the  original  samples;  this  is  much 
faster  than  other  types  of  interpolation  and  in  the  RAE  processor  the  intermediate 
samples  are  generated  via  the  FFT  technique  described  in  section  5.5  With  three 
intermediate  points  between  the  original  samples  the  ghosts  are  reduced  to  a 
level  of  about  -35  dB,  see  Fig  1 . 

16  CONCLUSION 

In  this  Report  an  attempt  has  been  made  to  set  down  the  theory  of  digital 
SAR  processing  for  orbital  SAR.  However,  the  reader  should  note  that  many  topics 
have  not  been  included.  For  example  there  are  several  other  azimuth  compression 
techniques.  In  this  connection  there  is  a  technique  which  at  RAE  is  called 
'coherent  subimaging'.  This  is  a  hybrid  time  domain-frequency  domain  technique 
in  which  a  large  number  of  prefilters,  say  16  or  32,  select  Doppler  bands  and 
an  image  is  generated  from  each  set  of  data.  The  resulting  images  are  not 
detected  but  are  added  coherently  to  produce  an  image  which  has  full  resolution. 
The  technique  has  several  variations,  the  point  being  that  because  the  data  is 
prefiltered  it  can  be  subsampled  and  all  the  benefits  of  data  reduction  obtained. 
Because  there  are  now  far  fewer  pulses  in  each  sub-aperture  the  time  domain 
approach  becomes  efficient  -  more  efficient  than  frequency  domain  processing  - 
see  the  comments  in  section  11  on  multilock  processing. 

Autofocussing  techniques  have  not  been  addressed,  nor  have  the  effects  of 
partial  coherence.  Another  problem  which  has  not  been  addressed  is  the  problem 
of  dynamic  range  in  SAR  processing,  that  is,  how  many  bits  does  one  need  at 
various  stages  of  the  processing  to  represent  the  data? 
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Many  of  these  tooics  are  under  active  investigation  at  RAE  and  this  list 
will  serve  to  demonstrate  that  the  study  of  methods  for  digital  SAR  processing  is 
a  problem  that  is  certainly  net  completely  solved;  much  work  is  still  needed 
in  many  areas,  for  example  on  achieving  higher  speed  processing,  on  selection 
of  optimum  radar  parameters  and  on  special  purpose  processing  to  enhance 
particular  image  features. 
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Fig  1  Image  of  point  target  (Goldstone  corner  reflector  from  SEASAT  rev.882) 
showing  the  difference  between  along  track  and  azimuth  directions 
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Figs  11-14 


Frequency 


Fig  1 1  Line  range  spectrum  caused  by  periodic  pulsing  of  the  radar 


Doppler  shift 


Frequency 


Fig  12  Line  split  due  to  Doppler  shift  with  cosine  phase  term 


Doppler  shift 
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Fig  13  Line  shift  due  to  Doppler  shift  with  exponential  phase  term 
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Fig  14  Spectral  lines  broadened  into  a  Doppler  continuum 
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Fig  15  Spectrum  of  interpolated  data 
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Figs  21-23 


Fig  21  Data  bananas 


Fig  22  Corresponding  spectrum  bananas 
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Fig  23  Piecewise  approximation  of  the  spectrum  banana 


Figs  24&25 


Fig  25  Piecewise  approximation  of  the  data  banana 
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Fig  27&28 
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Fig  27  Amplitude  modulation  of  the  data  due  to  piecewise  approximation 
in  the  time  domain 
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Fig  28  Amplitude  modulation  of  the  spectrum 
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